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Abstract 



A prescription is developed for matching general relativistic perturbations 
across singularities of the type encountered in the ekpyrotic and cyclic scenar- 
ios, i.e., a collision between orbifold planes. We show that there exists a gauge 
in which the evolution of perturbations is locally identical to that in a model 
space-time (compactified Milne mod Z%) where the matching of modes across 
the singularity can be treated using a prescription previously introduced by 
two of us. Using this approach, we show that long wavelength, scale- invariant, 
growing-mode perturbations in the incoming state pass through the collision 
and become scale-invariant growing-mode perturbations in the expanding hot 
big bang phase. 

PACS number(s): 11.25.-w,04.50.+h, 98.80.Cq,98.80.-k 



I. INTRODUCTION 



The big bang singularity is one of the most vexing puzzles in modern cosmology. Trac- 
ing time backwards, the field equations of general relativity break down in an apparently 
irretrievable manner some fourteen billion years ago when the density of matter and the 
curvature of space-time diverge. Cosmic inflation does not ameliorate this disaster, but 
rather tempts us to ignore it by just assuming that the universe somehow emerged from 
the singularity in an inflationary state, and that subsequent inflation washed out all of the 
details of the big bang and how inflation began. 

A more fundamental point of view is that the singularity is a manifestation of the break- 
down of general relativity at short distances, which needs to be properly dealt with in a more 
consistent cosmology. String theory and M-theory are important suggestions as to what a 
more fundamental theory might look like, improving on general relativity, for example, by 
providing consistent perturbative S-matrices that include graviton processes. If string the- 
ory is a consistent, unitary S-matrix theory, as it is believed to be, then it is reasonable 
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to expect that the cosmic singularity should be resolved within string theory, or a future 
development of it, in a satisfactory way. In particular, for every 'out' state there should 
be at least one 'in' state. The question arises: What could the 'in' state have been which 
produced the hot big bang? 

In recent papers, we have explored a concrete, detailed proposal for answering this deep 
question. In the ekpyrotic 1 and cyclic 2 Universe models, the origin of scale invariant density 
perturbations and the flatness, homogeneity and horizon puzzles of the standard cosmology 
are all explained without recourse to a burst of high energy primordial inflation. 1 ' 2 Instead, 
these puzzles are solved by physical processes occurring prior to the hot big bang, 1-4 in 
a highly economical way employing today's observed cosmological constant in an integral 
manner. However, key to the success of these new scenarios is a consistent passage through 
the big bang singularity. 

At first sight, passing safely through a big crunch/big bang transition seems impossi- 
ble because many physical quantities (density, curvature) diverge there. However, in the 
situation encountered in the ekpyrotic and cyclic brane world models, the situation is far 
less severe. 1 ' 2 When two boundary branes collide, even though this is the big bang sin- 
gularity in the conventional (Einstein frame) description, in the background solution the 
density of matter and the space-time curvature of the branes remain finite. Conservation 
of total energy and momentum across the collision may be consistently imposed 2 and, once 
the densities of radiation and matter generated on the branes at the collision are fixed (by 
microscopic physics), the outgoing state is uniquely determined. 

However, while the background geometry describing a boundary brane collision seems to 
be well behaved, it is still mathematically singular in the sense that one dimension disappears 
at one instant of time. The space-time ceases to be Hausdorff, 6 and since the dimensionality 
of the spatial slice is only three at this moment, it is not a good Cauchy surface. More 
worryingly, perturbations generally diverge as one approaches the singularity, as the result of 
the cosmological blue shift associated with the collapse of the extra dimension. Nevertheless, 
the situation is more manageable than it appears to be at first sight. In certain gauges, 
the metric perturbations only diverge logarithmically in time, 5 and the canonical momenta 
associated with the perturbations and certain other perturbation variables actually remain 
finite at the singularity. 

Around the brane collision, the space-time geometry may be modeled by a simpler space- 
time which we shall refer to as 'compactified Milne mod Z-i '. This is locally flat away from 
the singularity, and may be embedded within Minkowski space-time as shown in Figure 1. 
The model space-time may be thought of as describing the collision of two tensionless Z 2 
branes separated by a flat bulk. In a study of free fields on this space-time, two of us earlier 
showed 8 that the construction of a unitary map between incoming and outgoing states is 
not only possible but essentially unique. As we review in Section II, the basic idea is to 
employ normal propagation of free fields on the Minkowski covering space-time. This rule 
was shown 8 to satisfy many desirable properties. For example it defines a vacuum two-point 
function of Hadamard form which is also time reversal invariant. And in this idealized 
situation with no interactions, there turns out to be no particle production associated with 
passage through the singularity. Some first steps were taken towards studying interactions 
and these were shown to lead to finite answers provided the coupling constant vanishes 
sufficiently rapidly near the collision event. 
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FIGURES 




FIG. 1. Locally, the collision of two branes may be embedded in Minkowski space-time. The 
usual Minkowski space-time coordinates T and Y are are expressed as T = t coshy and Y = t sinhy, 
where the Lorentz-invariant coordinate t is constant on the dashed lines. The collision event is 
constructed in two steps. First the y coordinate is compactified by identifying y with y + 2yo, to 
produce the double-conical space-time shown at the right. Second, the circular sections of these 
cones are orbifolded by the Z2 symmetry y — > 2yo — y. The two fixed points of the Z2 symmetry 
are two tensionless branes moving at a relative speed of tanlryo, which collide and pass through 
one another at t = 0. 

The purpose of the present paper is to extend these ideas to a study of full general 
relativistic perturbations in space-times possessing singularities of the type shown in Figure 
1. The usual definition of a space-time manifold is that it is a metric space which appears 
locally flat. This means that in the neighborhood of any point P it should always be possible 
to choose a coordinate system in which (a) the metric at P is the Minkowski metric, and 
(b) the first derivatives of the metric with respect to each coordinate vanish at P. The 
inclusion of singular points of the type shown in Figure 1 requires an extension of these 
rules. In particular, the usual notion of general coordinate invariance becomes more subtle. 
A description of the incoming and outgoing space-times, away from the singularity, should be 
completely independent of coordinates since only the intrinsic geometry matters. However, 
connecting the two halves of the space-time across the singularity requires a correspondence 
between the 'incoming' and 'outgoing' coordinate systems. What this means in practice is 
that after solving for the metric and brane perturbations using general relativity in the upper 
and lower halves (which may be done in any gauge), one needs to choose a set of coordinates, 
or gauge, common to both halves within which the matching is to be performed. 
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FIG. 2. The definition of a space-time manifold is that when viewed 'up close' (left figure), it 
should appear to be locally flat. We define singular space-times of the type we are interested in 
here as space-times for which there exists a single coordinate system covering the neighborhood 
of the singularity in both the incoming and outgoing space-times, within which the collision event 
appears locally identical to the idealized situation of tensionless Z2 branes colliding in Minkowski 
space-time (right figure). 

Our proposal for extending general relativity to this type of singularity is illustrated in 
Figure 2. The idea is to insist that the the upper and lower halves be connected within a 
smooth set of embedding coordinates within which the geometry appears locally identical to 
that describing the model space-time consisting of the collision of two tensionless branes i.e. 
compactified Milne mod Zi- This set of embedding coordinates, locally unique up to Lorentz 
transformations, connects the contracting and expanding phases on either side of the bounce. 
The fact that fields may propagate across the singularity in the model space-time shown in 
Figure 1 and, at the same time, unitarity and all the other desirable physical properties 
of massless fields propagating in ordinary Minkowski space-time can be maintained 8 makes 
this minimal extension of general relativity that we propose both reasonable and physically 
sensible. 

In close analogy with the definition of a space-time manifold, we shall define 'locally' by 
insisting that the first two terms in a series expansion of the metric perturbations (specifi- 
cally the constant and logarithmic terms) behave precisely as free gravitational waves would 
in a compactified Milne mod Z2 space-time. The main work of the paper will be to demon- 
strate that this condition may be precisely formulated, at least for the lowest energy modes, 
and that it completely fixes the power series expansion in the Lorentz-invariant distance 
t = \/T 2 — Y 2 from the singularity. Within the coordinate systems so constructed for the 
incoming and outgoing space-times, we find a unique rule for matching gravitational per- 
turbations, in a manner entirely analogous to the matching of free scalar fields in the model 
space-time, as discussed in Ref. 8. 

The matching procedure we propose is more subtle than that usually adopted in general 
relativity. In situations where the matter stresses change suddenly on some physically pre- 
scribed space-like surface (for example in a phase transition), it is normally only necessary 
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to match the spatial three-metric and its normal time derivative, without worrying about 
the detailed behavior of the solutions of the field equations. In our case, the metric pertur- 
bations diverge at the singularity. One might attempt to cut the divergence off by pasting 
the incoming to the outgoing space-time together on some arbitrary surface slightly away 
from the singularity, but it is not known how to do this in a coordinate-invariant manner 
inevitably leading to ambiguous, and usually cutoff-dependent answers. In contrast, our pro- 
cedure for massless fields including gravitational waves on compactified Milne mod Z2 can 
be formulated in terms of analytic continuation, which is automatically coordinate- invariant, 
or in terms of a real continuation in an embedding Minkowski space with asymptotically flat 
boundary conditions, also a coordinate-invariant prescription. Both methods produce the 
same cutoff independent result. Notice also that both involve global aspects of the space- 
time, and cannot be stated as a purely local matching rule. This seems to be the inevitable 
price one has to pay for evolving through a singularity where a Cauchy surface does not 
exist. 

We have in mind of course, an application of this proposal to the types of cosmological 
singularities encountered in ekpyrotic and cyclic models in which two boundary branes collide 
as shown in Figure 1. In particular we wish to track scale- invariant perturbations developed 
via the ekpyrotic mechanism 1 ' 5 in the incoming state across the singularity and into the 
outgoing hot big bang phase. The conclusion of our work is that with the prescription 
adopted here, scale-invariant, growing mode perturbations produced during the pre-big bang 
phase 1 ' 2 ' 4 pass through the bounce and become scale-invariant growing mode perturbations 
in the late Universe. 

Let us briefly comment on the relation of this paper to previous studies by ourselves 
and others. Our first attempt 5 at matching perturbations across the transition was based 
entirely on the study of the four-dimensional effective theory. As we shall see, this is not 
sufficient to describe the bounce, which is really five- dimensional. Nevertheless, in that 
work we observed that certain perturbation variables, such as the comoving energy density 
perturbation e m were finite at the singularity and could be matched across it. The present 
(and far more sophisticated) approach confirms this element of the procedure. The problem 
is that two matching conditions are needed in the four-dimensional effective theory and the 
first derivative of e m turns out not to be independent of e m itself because the differential 
equation is singular at t — 0. This leads to an ambiguity in the second matching condition. 
Based on simplicity, we proposed matching the second derivative and obtained an outgoing 
scale- invariant spectrum. However, we did not have any real physical justification for this 
choice. 

There were criticisms and alternative proposals for matching conditions, 9 including the 
idea that one should match the curvature perturbation on comoving (or constant density) 
slices, 10-12 a procedure which is often useful in the context of nonsingular, expanding four- 
dimensional cosmology. In our setting, the comoving curvature perturbation is logarithmi- 
cally divergent at the singularity, 5 ' 17 but if one disregards this and proceeds to match its 
long wavelength, constant component, this proposal results in the growing, scale-invariant 
perturbations present in the pre-big bang phase being matched to a pure decaying mode in 
the outgoing state. 11 ' 12 The result is a complete absence of long wavelength density pertur- 
bations in the big bang phase. It was subsequently pointed out, however, that this null result 
is atypical in the sense that, for most choices of matching surfaces, scale invariant growing 
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perturbations coming in would match to scale invariant growing perturbations coming out. 16 
Some of the alternative proposals are designed specifically for four-dimensional theories 
in which the bounce from contraction to expansion occurs at a non-zero value of the scale 
factor 13 ' 14 (see also Ref. 15). This is accomplished by arranging for the equation of state w 
to violate the null energy condition near the bounce, i.e., w < —1. We emphasize that the 
ekyprotic and cyclic scenarios and the considerations here do not fall into this category. The 
bounce in Figure 1 corresponds to zero scale factor in the four-dimensional effective theory 
and the four- dimensional effective equation of state parameter is strictly positive before the 
collision. 

Our approach is to choose a class of gauges in which the geometry around the collision 
event appears locally identical to that describing linearized perturbations around the model 
space-time, compactified Milne mod Z2. Then we match the perturbations according to the 
procedure of Ref. 8 for that space-time. An important feature of our choice of coordinates 
is that the collision event is simultaneous in Milne time and occurs at the background value 
t — both for the incoming and outgoing state. That is, the limits t — > 0~ in the incoming 
state and t — > + in the outgoing state correspond to the same physical space-time surface. 

In the course of our analysis we shall uncover the problem with matching the curvature 
perturbation on comoving (or constant energy density) slices in the four-dimensional effective 
theory, (4, across the bounce. We shall show that £4 is indeed conserved on long wavelengths 
both before and after the bounce, and that furthermore on long wavelengths it is equal to 
the comoving curvature perturbations on the branes (±. Why then are these variables not 
conserved across the bounce? The reason, detailed in Section V.E, is that the brane collision 
event is not simultaneous in the comoving or constant energy density time slicing. This is a 
disaster in terms of matching. In this coordinate system, the t — > + and t — > 0~ space-like 
surfaces do not physically coincide and therefore perturbations should certainly not match 
across them. We find that the collision event is displaced from the t = + and t = 0~ surfaces 
in these slicings by a scale-invariant time delay, within which all the information regarding 
the growing mode perturbation is contained. A determination of the collision-synchronous 
time slices is only possible within the full five-dimensional theory, and our final result for the 
spectrum of growing mode perturbations involves five-dimensional parameters which cannot 
be re-expressed in purely four-dimensional terms. 

Distinct but closely related are problems raised in recent attempts to directly study string 
theory on compactified Milne space-times analogous to that shown in Figure 1. 20 > 21 Since 
these types of background are locally flat, one can solve 18 the tree level field equations of 
string theory to all orders in a', away from the singularity. It is then tempting to calculate 
string scattering processes using a Lorentzian generalization of standard orbifold techniques 
to this time-dependent case. Calculations have been performed in analogous backgrounds, 
for example, the null orbifold and 'null-brane' backgrounds 19,20 possessing some remaining 
super symmetry. The result is that tree level scattering amplitudes develop infrared diver- 
gences which have been attributed to the back-reaction of the geometry near the singularity. 

It is unclear what the physical significance of these results are yet. The breakdown of 
string perturbation theory seems to indicate that nonlinear effects must be taken into ac- 
count. But such nonlinear effects are not necessarily disastrous for cosmology. For example, 
since the collision takes place on a very short timescale, one plausible possibility is that non- 
linearities result in the production of microscopic black holes at the collision. This would 
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be consistent with the conclusion that perturbative string theory breaks down, but it would 
be unimportant for cosmology. The black holes would radiate and decay rapidly after the 
bounce without having a significant effect on the long wavelength perturbations that are 
relevant cosmologically. 

The classical theory may provide some insight. For example, consider classical general 
relativity with a scalar field. As the universe contracts towards a big crunch singularity, 
the gradients of the energy density diverge and one might be tempted to argue that the 
homogeneous Friedmann- Robert son- Walker (FRW) equations become invalid, however, this 
conclusion is believed to be wrong. Instead, the behavior of the metric and fields becomes 
ultralocal. 46 Spatial derivatives become less important as the the universe contracts and, at 
each point in space, the geometry follows a homogeneous evolution. This occurs because, 
although the gradient terms grow, the homogeneous terms grow faster. A description of 
this subtle situation may well be difficult using string perturbation theory, which relies 
for example upon the existence of a globally good gauge. However, as we shall explain 
in the conclusions, there is a simple classical picture of where the nonlinearities lead to. 
And within this picture, we see that the nonlinear corrections would hardly alter our final 
matching result. 

We should also note that the string theoretic calculations have only so far been possible 
in certain special models for which the technical tools needed are available. In particular, 
they have all been done in the context of ten dimensional string theory at fixed coupling, 
using Lorentzian orbifolding, with one of the nine string theory spatial dimensions shrinking 
away and reappearing. However, this setup is quite different from the case proposed for the 
ekpyrotic model, where the tenth spatial dimension (of eleven dimensional supergravity), 
separating the two boundary branes, was supposed to collapse and reappear. The eleven 
dimensional theory reduces, at fixed, small brane separation, to string theory 23 at weak 
coupling. But in the time- dependent situation we are interested in, the coupling would 
actually vanish as the branes meet. This situation is qualitatively different from the examples 
which have been studied so far. In particular, the infinities encountered in Refs. 20 are 
proportional to the string coupling. But in the ekpyrotic model the coupling vanishes at the 
singularity. 

Progress in the investigation of such singularities within string theory 31 continues to 
be an active field. 22 ' 24-27 Analytic continuation methods related to those we employed for 
field theory 8 have been applied to constructing string theory on similar backgrounds 28 with 
less pessimistic conclusions than the above cited works. 29 Other approaches and methods 
have also been developed. 30 We have continued to develop a simpler field theoretic approach, 
because it is considerably more manageable and may yield helpful physical insight. We hope 
that further developments of string theory can be used to check and develop the approach 
presented here. 

The remainder of the paper builds in stages towards a full calculation of the propagation 
of cosmological perturbations through a bounce of the ekpyrotic/cyclic type: 

• We first consider the propagation of scalar fields in a fixed background corresponding 
to two tensionless Z 2 branes colliding in a flat bulk as discussed in Ref. 8. [Section II] 

• We next consider linearized gravitational perturbations of the same model space-time. 
[Section IV] 
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• Finally, we consider the full-blown calculation of cosmological perturbations for two 
colliding branes with tension and a warped bulk. This calculation leads to our central 
result for the amplitude of the scale-invariant perturbations propagating across the 
singularity into the hot big bang phase. [Section V] 

Various tools are developed along the way. Section III develops the moduli space approx- 
imation for two colliding branes in a negative cosmological constant bulk 35 which we shall 
study as our canonical example. We extend this formalism, showing for example that it is 
exact for empty branes at arbitrary speed and curvature. In Appendix 1 we show that the 
four-dimensional effective theory consistently predicts the projected Weyl tensor contribu- 
tion to the effective Einstein equations on the branes, and is in agreement with the recently 
developed 'covariant curvature' approach 43 as well as earlier metric based approaches. 41 ' 42 
We also match the parameters of four-dimensional effective theory for the homogeneous flat 
background solution to the parameters of the five-dimensional theory. Appendix 2 discusses 
the gauge invariant variables for the five-dimensional theory and how the position of the 
branes depends on the choice of gauge. Appendix 3 works out the detailed background ge- 
ometry near the bounce in a coordinate system convenient for the perturbation calculations. 
Appendix 4 concerns the choice of gauge required to have the brane collision simultaneous 
at all values of the noncompact coordinates x 



II. PROPAGATION OF SCALAR FIELDS IN A COLLISION OF TENSIONLESS 

BRANES 

The idealized space-time we shall use as a model for the singularity is just Minkowski 
space-time subject to two identifications. 7 Expressing the usual Minkowski coordinates as 
T — t coshy and Y — t sinhy, the line element is 

ds 2 = -dT 2 + dY 2 + dx 2 = -dt 2 + t 2 dy 2 + dx 2 . (1) 

The incoming and outgoing regions, respectively t < and t > 0, are the two halves of Milne 
space-time Ai x R 3 . We now compactify the y coordinate by identifying under boosts, which 
correspond to translations in y, y — > y + 2y . We refer to the resulting space as compactified 
Milne space-time, or M c x R 3 . Finally we introduce two tensionless Z 2 branes by identifying 
fields under reflection across the circle, y — > 2y — y giving the orbifolded space Ai c jZ 2 x R 3 , 
or compactified Milne mod Z 2 . The branes are separated by a coordinate distance Ay = y 
which is the rapidity associated with their relative speed. Later in the paper it will be 
convenient to choose a Lorentz frame in which the branes are located at equal and opposite 
values of y — ±y /2. Note that any field which is even under the Z 2 must obey Neumann 
boundary conditions d y ip = on the two branes. 

The problem of propagating a free quantum field through a big crunch/big bang singu- 
larity of the type shown in Figure 1 was considered in Ref. 8. The equation of motion for a 
scalar field on the background (1) is 

1 k 2 

<P+-<P + -fi<P + k 2 <p = 0, (2) 
where k y is the momentum in the y direction and k that in the uncompactified x directions. 
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In this paper, our main interest is in the lowest excitations corresponding to the modes of 
the four-dimensional effective theory. In this compactified Milne setup these modes are the 
y-independent fields, trivially satisfying Neumann boundary conditions on the branes and 
periodicity in y. For these modes, equation (2) is just Bessel's equation with index v — 0. 
The two linearly independent solutions are Jo(kt) and N (kt), behaving for small positive t 
as 

J Q (kt) ~ 1 + ..., N (kt) ~ — (ln(/ct) + 7 — ln2) + . . . , (3) 

7T 

where 7 is Euler's constant 0.577. . .. The positive (respectively negative) frequency outgoing 
modes ip^ {i>^) are those which tend to the adiabatic positive (negative) frequency solu- 
tions as t — > 00. They are proportional to the Hankel function = J — iN (respectively 
Hq 1 ^ = J + iN ), and converge rapidly to zero in the lower (upper) half complex t-plane. If 
we split the quantum field tp(t, x) into its positive and negative frequency parts, they are well 
defined respectively in the lower and upper half complex t-plane. The unique analytic con- 
tinuation from negative to positive values of t is then to continue the positive frequency part 
below and the negative frequency above the singularity ait — 0. Continuing the expressions 
(3) around a small semicircle below t — one infers the relation H^\kt) = —H^\—kt) 
giving the positive frequency mode function at negative values of t. We can translate this 
into a matching rule for the field by writing Lp = £)a^ + ) + h.c, with a arbitrary and 
complex. The asymptotic behavior of the field (f is then found to be 

V ~ Qin + PirMk\t\ >0~, ~ Q out + P out lnk\t\ t -> + , (4) 

and the above continuation implies that 

Qout = —Qin + 2(7 — hl2)Pj n , P^t = Pi n . (5) 

The canonical momentum of the field \t\(p is actually proportional to sign(t)P. Hence, the 
field momentum reverses at t — with this matching rule. Note, however, that the constant 
term Q is not preserved across t — 0. Hence this matching rule is not simply time reversal 
at t — 0, and there is an arrow of time across t — 0. 
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FIG. 3. Continuation of left and right moving modes. A free field propagating in the lower 
quadrant may be decomposed into left and right movers as it approaches the past light cone of 
the origin T = Y = 0. The left movers are regular across Y = T < and may be continued 
into the left quadrant Y < 0, \T\ < The right movers are regular across the right segment 
Y = —T > and may be continued into the right quadrant Y > 0,\T\ < \Y\. If we impose 
vanishing boundary conditions at large Lorentz-invariant separation from the origin in the left and 
right quadrants, then once we know the left mover in the left quadrant, the right mover on the null 
segment Y = —T < is uniquely determined, and similarly the left mover on Y = T > 0. One 
thereby obtains a unique matching rule from the incoming, lower quadrant to the outgoing, upper 
one. 

There is another way of looking at this rule which is illustrated in Figure 3. Take a 
field configuration on one copy of the incoming wedge and repeatedly reflect it through 
the boundary branes to fill out the lower quadrant. The resulting configuration obeys the 
field equation (even with nonlinear interactions), as long as the equation is Z 2 invariant. 
The solutions to the field equation then naturally split into left and right movers as one 
approaches the light cone. The left movers are regular on Y = —T and the right movers 
on Y = T. Each can therefore be uniquely matched across the appropriate segments of the 
past and future light cone of the singularity (Figure 3). 

In this way, incoming data in the lower quadrant uniquely determines the left moving 
modes entering the left quadrant and the right moving modes entering the right quadrant. 
The solutions in the left and right quadrants may be fully specified by choosing boundary 
conditions. It is natural to demand that the fields vanish at space-like infinity. Once the 
solution in the left and right quadrants is determined then the left movers from the right 
quadrant and the right movers from the left quadrant may be uniquely matched to the 
left and right movers in the upper quadrant, completely determining the solution in the 
outgoing state. Again, in the context of our model spacetime compactified Milne mod Z 2 , 
this prescription yields exactly the same matching rule (5). The advantage of this derivation 
is that it gives the clearest explanation for the sign change in the constant contribution Q, 
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between the 'in' and 'out' states. This is just due to our having imposed a 'reflecting' 
boundary condition at space-like infinity. Since in passing from the lower to the upper 
quadrant, one such reflection is involved, a relative minus sign is acquired. And as we shall 
explain in Section IV, precisely the same matching rule may be applied for gravitational 
perturbations on compactified Milne mod Z 2 - In this case one can see that the condition of 
asymptotic flatness imposed in the two unphysical quadrants is actually coordinate invariant. 

In the case of cosmological interest where the branes have tension and the bulk is warped, 
the sign change of Q in (5) is still guaranteed provided two reasonable conditions are fulfilled. 
Assume that the low energy modes in the space-like regions (which are just the analytic 
continuation of the corresponding modes in the lower quadrant, obtained by setting t = is 
and y = p — in/2, where T = ssinhp and Y = scoshp), depend only on s as s — > 
(i.e. behave as the Kaluza-Klein zero modes). Second, assume that the mode selected by 
the imposed boundary condition at spacelike infinity behaves, near s = 0, as D + ln(fc|s|) 
with D a model-dependent constant. This is the generic behavior - for compactified Milne 
mod Z 2 we have D = 7 — ln2. Then it is straightforward to show by explicit calculation 
that matching the left/right movers across the light cone from the lower quadrant into 
the left/right quadrants and then into the upper quadrant, one obtains P ou t = Pi n and 
Qout = —Qin + 2-DPj„. Hence we see the sign change of Q is universal but the coefficient D 
is not. 

It is important to emphasize that all of these arguments for the matching rule (5) involve 
the detailed global structure of the embedding space-time. In particular the 7 — In 2 term in 
(5) is peculiar to the Minkowski embedding spacetime appropriate for compactified Milne 
mod Z 2 - If the embedding space-time is warped, the corresponding constant would be altered 
to some constant D as explained above. Fortunately it shall turn out that for the case we 
are interested in, P in <C Q in at long wavelengths and hence we are insensitive to the value 
of D. The correspondence Q out ~ —Qin is however universal as argued above and therefore 
reliable even in the warped case. It turns out that this sign change is crucial in allowing 
scale invariant growing perturbations to propagate across the singularity, in the absence of 
radiation. Furthermore, the sign change is interesting and important in the nonlinear theory, 
as we explain in the conclusions. 

III. THE 4D EFFECTIVE THEORY 

In subsequent sections we shall extend the matching rule just discussed for free scalar 
fields to full general relativistic perturbations. There are two major complications. The first 
is the gauge invariance of general relativity which, as explained above, is unusually subtle 
for singular space-times such as we are dealing with. The second is that the bulk space-time 
is not globally Minkowski space-time but is warped and has non-negligible y-dependence. 
Of course, this is related via Israel matching (see e.g. Ref. 34) to the fact that the brane 
tensions are nonzero. 

We want to solve the linearized Einstein field equations for five-dimensional gravity 
coupled to a pair of colliding orbifold (Z 2 ) branes. For the cosmological applications, we 
need to follow the system from times well before the brane collision, when the scale-invariant 
perturbations were generated, through the collision and into the far future. In general this 
would involve solving a system of coupled partial differential equations in y and t for the 
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bulk gravitational fields with mixed boundary conditions following from the Israel matching 
conditions on the branes, and would be well beyond an analytic treatment. 

However, there is a powerful tool we can call upon which makes the task surprisingly 
tractable: the moduli space approximation. 



On general grounds one expects the long wavelength, low energy modes of the system 
to be described by a four-dimensional effective theory, and we are only interested in low 
energy incoming states which are well described by this theory. We shall show that the four- 
dimensional effective theory may be consistently used to predict the brane geometries all the 
way to collision, thereby providing boundary data for the bulk five-dimensional equations 
which we solve as an expansion in t about the collision event. After the collision, the 
four- dimensional effective theory plays an equally important role, enabling us to track the 
behavior of perturbations into the far future of the collision event (Figure 4). The technique 
we describe forms the basis for our analysis of the singularity described in later sections, 
but it is also of considerable generality and use in its own right, since almost all of the late 
Universe phenomenology of brane worlds can be most efficiently described using the effective 
theory alone. 



FIG. 4. The worldlines of the positive and negative tension branes are plotted for some fixed 
value of the uncompactified coordinates x. The four-dimensional effective theory is used to predict 
the intrinsic geometries of the positive and negative tension branes, i.e. their space-time metrics 
<7+, and g~ u , according to equation (10). The four-dimensional effective theory is used to describe 
the incoming and outgoing perturbed branes far to the past or future of the collision event. The 
brane metrics also provide boundary data for the five-dimensional bulk metric which we solve for 
as a power series expansion in time about the collision event. 

In this paper we concentrate on the simplest two-brane world model consisting of one 
positive and one negative tension brane bounding a bulk with a negative cosmological con- 
stant A = — 6M|/L 2 where L is the AdS radius and M 5 the five-dimensional Planck mass. 
If the brane tensions a± are fine tuned to the special values ±6M|/L, the system allows 



A. The Moduli Space Approximation 
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a two-parameter family of static solutions in which the scale factor on each brane is a free 
parameter, or modulus. The idea of the moduli space approach is that such parameters are 
promoted to space-time dependent fields within the four-dimensional effective theory. In 
passing, we note that many of the methods we use in this paper should in principle extend 
to more complicated theories such as Horava-Witten theory, in which the family of static 
solutions exists without the need for fine tuning of the brane tensions. 

In Khoury et al., 1 the effective action for the moduli in this system was computed in the 
low velocity approximation, and shown to be equivalent to Einstein gravity plus a scalar 
field which couples non-minimally to the matter on each brane (see also Ref. 36). The 
derivation given here, while more specific to the simplest brane models, is both simpler 
and more powerful. It shows that the same effective action actually has a broader range of 
validity than originally anticipated, turning out to be exact for empty brane configurations 
with cosmological symmetry, for arbitrary spatial curvature and velocity (or expansion rate). 
When matter is present, the effective theory is a good approximation as long as the density of 
matter is small compared to the brane tension. The fact that the four-dimensional effective 
theory is so accurate is likely to be a special feature associated with the lack of bulk degrees 
of freedom in the simplest brane world model we are focusing on: for configurations with 
cosmological symmetry, a generalized Birkhoff theorem 37 ' 40 holds which guarantees that no 
radiation is emitted into the bulk. 

Consider a positive or negative tension brane with cosmological symmetry but which 
moves through the five- dimensional bulk. The motion through the warped bulk induces 
expansion or contraction of the scale factor on the brane. As shown in Ref. 34, the scale 
factor on the brane obeys a 'modified Friedmann' equation, 

2 1 p\ K c 

where p± is the density (not including the tension) of matter or radiation confined to the 
brane, b± is the brane scale factor, and H± is the induced Hubble constant on the positive 
(negative) tension brane. We work in units such that the coefficient of the Ricci scalar in 

M 3 

the five-dimensional Einstein action is The last term is the 'dark radiation' term, where 
the constant C is related to the mass of the black hole in the Schwarzchild-AdS solution 
discussed in Appendix 3. 

We shall show that the solutions to these equations are precisely reproduced by a four- 
dimensional effective theory, with the only approximation necessary being that the density 
of matter or radiation confined to the branes, p± be much smaller than the magnitudes of 
the brane tensions, so that the p\ terms in (6) are negligible. For the particular concerns in 
this paper, namely the accurate calculation of the long wavelength curvature perturbation 
on the branes, it is reassuring that the four-dimensional effective theory description is such 
a well-controlled approximation, even at large brane velocities, in the long wavelength limit. 

Choosing conformal time on each brane, and neglecting the p 2 terms equations (6) become 
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where prime denotes conformal time derivative. The corresponding acceleration equations 
for and b'!_, from which C disappears, are derived by differentiating equations (7) and using 
d(pb 4 ) = b 3 (p — 3P) db, with P being the pressure of matter or radiation on the branes. We 
now show that these two equations can be derived from a single action provided we equate 
the conformal times on each brane. Consider the action 



/ 



dtNd 3 x 



-ZMlL(N- 2 b'l - Kb%) - p+b% + 3L(N~ 2 b'f - Kb 2 + ) - p.b 



(8) 



where N is a lapse function introduced to make the action time reparameterization invariant. 
Varying with respect to b± and then setting N — 1 gives the correct acceleration equations 
for 6" and b"_ following from (7). These equations are equivalent to (7) up to two integration 
constants. The constraint equation, following from varying with respect to iV and then set- 
ting N — 1, is just the difference of the two equations (7) and ensures that one combination 
of the integration constants is correct. The constant C is then seen to be just the remaining 
constant of integration of the resulting system of equations and can in effect be determined 
by the solutions of equations of motion following from the action (8). 38 

Having shown that the modified Friedmann equations (with the neglect of p 2 terms) 
follow from an action in which C does not appear, we are now able to change variables to 
those in which the system appears as conventional Einstein gravity coupled to a scalar field 
plus matter. We rewrite the action (8) in terms of a four-dimensional effective scale factor 
a and a scalar field 0, defined by b + = a cosh(0/v^6), b- = — asinh(0/v^6). Clearly, a and 
transform as a scale factor and as a scalar field under rescalings of the spatial coordinates 
x. To interpret more physically, note that for static branes the bulk space-time is perfect 
Anti-de Sitter space with line element dY 2 + e 2Y ^ L (— dt 2 + dx 2 ). The separation between the 
branes is given by d — L\n(a + / a ) = Lin (— coth(0/-\/6)), so d tends from zero to infinity 
as tends from minus infinity to zero. 

In terms of a and 0, the action (8) becomes 



dtd x 



-3M 3 L(a 2 - Ka 2 ) + -a 



2/2 



(9) 



which is recognized as the action for Einstein gravity with line element a 2 (t)(—dt 2 + 
r y i jdx l dx : >), being the canonical metric on H 3 , S 3 or E 3 with curvature K, and a min- 
imally coupled scalar field 0. The matter action S m is conventional, except that the scale 
factor appearing is not the Einstein-frame scale factor but instead b + = a cos h(0/ v / 6) and 
6_ = — a sinh(0/v^6) on the positive and negative tension branes respectively. 

Now we wish to make use of two very powerful principles. The first is the assertion that 
even in the absence of symmetry, the low energy modes of the five- dimensional theory should 
be describable with a four-dimensional effective action. The second is that since the original 
theory was coordinate invariant, the four dimensional effective action must be coordinate 
invariant too. Since the five-dimensional theory is local and causal, it is reasonable to 
expect these properties in the four- dimensional theory. If furthermore the relation between 
the four-dimensional induced metrics on the branes and the four-dimensional fields (i.e. the 
four-dimensional effective metric and the scalar field 0) is local (as one expects for the long 
wavelength, low energy modes we are interested in), then covariance plus agreement with 
the above results forces the relation to be 
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g+ = (cosh^/VetfgX v = (-mnhtf/VE)) 2 g% . (10) 

When we couple matter to the brane metrics, these expressions should enter the action 
for matter confined to the positive and negative tension branes respectively. Likewise we 
can from (9) and covariance immediately infer the effective action for the four- dimensional 
theory: 

r I M 2 1 \ 

S = j d\^-g \-fR - -(«9,0) 2 J + S-[g-} + S+[g + ], (11) 

where we have defined the effective four- dimensional Planck mass Mf = (87^4) _1 = MfL. 

B. Branes with non-zero matter density 

For most of this paper we shall only study the specially simple case of radiation on the 
branes (which are 3+1 dimensional). The matter action is then independent of as a result 
of the conformal invariance of radiation in 3+1 dimensions, and this will greatly simplify 
our analysis. But as an aside let us for a moment consider nonrelativistic matter on the 
branes. Then there is a non-minimal coupling with 0, leading to a source term in the scalar 
field equation: 

□0 = -i(cosh( A)4) _ 1 ( Sinh (-|)^T-, (12) 

where primes denote derivatives and the are the traces of the stress tensors for matter 
on the two branes contracted with respect to the relevant brane metric. It is interesting 
to see how these results compare with what is known about brane world gravity from prior 
studies. 39 For perfect fluids, the effective matter Lagrangian 32 reads — / d 4 x\/—g ± p±. Hence 
matter on the branes couples to the four-dimensional (Einstein frame) effective theory in 

the combination = cosh(0/v^6) p+ + sinh(0/v^6) 4 p-- As the inter-brane distance grows, 
the field tends to zero. Since the cosh tends to unity, we see that a matter source on the 
positive tension brane with physical density p + contributes the same amount to the density 
seen by Einstein gravity in the four- dimensional effective theory. Furthermore, from (12), 
the coupling of such matter to the dilaton vanishes as 0. Hence the dilaton decouples and 
ordinary Einstein gravity is reproduced in this limit. Matter on the negative tension brane 
behaves very differently. If its density as seen by Einstein gravity in the four- dimensional 
effective theory is p 4 , then its physical density on the brane is much larger, p_ ~ 0~ 4 P4, and 
from (12) it sources the dilaton field as _1 p4- Hence at small the source for the dilaton 
diverges and Einstein gravity is never reproduced. 

The derivation we have just given of the four- dimensional effective action starting from 
the modified Friedmann equations is in the present context both simpler and more powerful 
than previous derivations. It shows that the induced geometries on the branes are correctly 
predicted for branes with cosmological symmetry, for arbitrary curvature and speed of the 
branes provided only that that the p 2 matter terms are negligible. For these cosmological 
backgrounds, the four-dimensional effective theory accurately describes the brane collision 
even though from the Einstein frame point of view such a collision is highly singular in the 
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sense that the 4d effective scale factor a tends to zero, the Riemann tends to minus infinity 
in finite time. Nevertheless, the brane geometries and densities described by g^ v and p±, are 
finite and well behaved at all times. 

One surprising point about the map from five- dimensions to four is that the effective 
theory with a scalar field sourced by the combined energy density p 4 = cosh(0/v^6) p+ + 
sinh(0/\/6) 4 p- manages to correctly predict the solutions to the Friedmann equations on 
each brane even though these are separately sourced by p + and p_. This is possible because 
of the integration constants. In the four-dimensional effective theory the basic equations 
can be taken to be the Friedmann equation (a' 2 = . . .) which has one integration constant 
and the scalar field equation ((a 2 0')' = . . .) which has two. So there is a three-parameter 
set of solutions, although one of these is not physical as it is just a rescaling of a. On the 
other hand the two brane Friedmann equations have two integration constants along with 
the additional constant C which is the dark radiation term. Consequently we have a precise 
match between the integration constants showing that there is a one to one map between 
the solutions of the two sets of equations. In performing an explicit check we find that the 
missing information on how much matter is contained on each brane is contained in the 
integration constants for the dilaton equation. 

C. Relation between 4d effective theory and 5d brane parameters 

The five-dimensional background we seek to describe consists of two parallel, flat Z 2 - 
symmetric three-branes bounding a bulk with a negative cosmo logical constant. In the 
incoming state, as they head towards a collision, the branes are assumed to be empty. In 
the ekpyrotic scenario, it is assumed that the brane collision event fills them with radiation. 
In this section we shall see how to describe this background setup in terms of the four- 
dimensional effective theory, and in particular we shall determine precise relations between 
the parameters of the four and five-dimensional theories. The two brane geometries are 
determined according to the formulae (10), and the background solution relevant post- 
collision is assumed to consist of two flat, parallel branes with radiation densities p±. The 
corresponding four- dimensional effective theory has radiation density p r , and a massless 
scalar field with kinetic energy density p^. It is convenient to work in units where the four- 
dimensional reduced Planck mass M 4 = (8irG)~z is unity. The four-dimensional Friedmann 
equation in conformal time then reads 

o* = h pr a 4 + P0 a 4 ) = 4A 4 (r 4 + ^), (13) 

o (X 

where we have defined the constants A 4 and r 4 , and used the fact that the massless scalar 
kinetic energy p^ oc a~ 6 . The reason for this choice of constants will become clear momen- 
tarily. 

The solution to (13) and the massless scalar field equation (a 2 0')' = is: 

« J = ^(l + r 1 r), ^=/| ta ((Tf^))' (14) 

From these solutions, we reconstruct the scale factors on the branes according to (10), 
obtaining: 
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b± = 1 ± A 4 r + r 4 r, (15) 

so we see that with the choice of normalization for the scale factor a made in (13), the brane 
scale factors are unity at collision. For comparison, in Ref. 2 we parameterized the radiation 
density appearing in the four dimensional effective theory using the Hubble constant H r 

at equal density of the radiation and scalar kinetic energy, H r = (2r 4 )2 / A\. Also, the 
parameter H 5 used there to describe the contraction rate of the fifth dimension may be 
expressed, for r ± L 2 « 1 and slow velocities as 2^4 4 . 

We may now directly compare the predictions (15) with the exact five-dimensional solu- 
tion given in in equations (123) of Appendix 3, equating the terms linear in r to obtain 

A 4 = (1/L)(1 + L2(r+ 1 ~ r - ) ) tanh(y /2), 

n = L(r+ + r - } (16) 

4 12tanh(y /2)' 1 ' 

where y is the rapidity associated with the relative velocity of the branes at collision V = 
tanh(y ) and r± is the value of the radiation density p± on each brane at collision. These 
formulae are the exact expressions for the four-dimensional parameters in terms of the five- 
dimensional parameters neglecting contributions of order p 2 . In fact, at leading order in r 
they are better than this since to this order the four-dimensional prediction is exact. 

For later purposes it will also be useful to define the fractional density mismatch on the 
two branes as 

/ = r -^, (17) 



so that we have 



12 fr A 

r+ - r_ = — ^ tanh(y /2). (18) 



D. Four Dimensional Perturbation Equations 

In this section, we describe the perturbations of the brane-world system in terms of the 
four-dimensional effective theory. The only cases we consider in detail are where the branes 
are empty or carry radiation. The conformal invariance of radiation in four-dimensions 
greatly simplifies matters since the scalar field then has no direct coupling to the radiation 
and hence the latter evolves as a free fluid in the four-dimensional effective theory. We 
elaborate on the significance of this conformal invariance in section VI, part C. 

We shall now describe the scalar perturbations, in longitudinal (conformal Newtonian) 
gauge with a spatially flat background where the scale factor and the scalar field are given 
by (14). The perturbed line element is 

ds 2 = a 2 (r) (-(1 + 2$)dr 2 + (1 - 2tf) dx 2 ). (19) 

Since there are no anisotropic stresses in the linearized theory, we have $ = \l> (see e.g. 
Ref. 32). 
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A complete set of perturbation equations consists of the radiation fluid equations, the 
scalar field equation of motion and the Einstein momentum constraint: 



6' r = ~(k 2 v r -W) 

v' r = -5 r + $ 
{5<j))" + 2H{5<t>)' = -k 2 (5(j)) + 40'$ 

& + H$ = la 2 p r v r + ±<j>'(8<f>), 



(20) 



where primes denote r derivatives, 5 r is the fractional perturbation in the radiation density, 
v r is the scalar potential for its velocity i.e. v r = Vi> r , S(p is the perturbation in the scalar 
field, and from (14) we have the background quantities H = a' /a = (l + 2r 4 r)/(2r(l + r 4 r)), 



We are interested in solving these equations in the long wavelength limit, \kr\ <C 1. There 
are only two independent solutions to (20), namely a growing and a decaying mode, provided 
that we specify that the perturbations are adiabatic. Recall that the idea of adiabaticity in 
the cosmological context is that for long wavelength perturbations, there should be nothing 
in the state of the matter to locally distinguish one region of the Universe region from 
another. At each spatial location the evolution of the densities of all the different fluids 
(radiation, baryons, dark matter) should a single history in which each fluid evolves with 
the scale factor a according to dpi = — 3(pi + Pi) dlna = — 3pj(l + Wi) dlna where pi is its 
density, Pi is its pressure and Wi parameterizes the equation of state. Likewise the total 
density evolves as dp = —3{p + P) dlna = — 3p(l + w) da. Since the history is parameterized 
uniquely by the scale factor a, an adiabatic perturbation can be thought of as arising from 
a fluctuation 5 In a. Hence solving all the above equations for 5 In a, one finds 



for adiabatic perturbations. 

For the case at hand, the components of the background energy density in the four- 
dimensional effective theory are scalar kinetic energy, with — 1, and radiation, with 
w r — |. It follows that for adiabatic perturbations, at long wavelengths we must have 



In longitudinal gauge, the fractional energy density perturbation and the velocity potential 
perturbation in the scalar field (considered as a fluid with w = 1) are given by 



and > /|^ = 1/(t(1 + r 4 r)). 




(22) 




(23) 



From the equations (20) above (and using <p' oc a 2 ) it follows that 




(24) 
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Maintaining the adiabaticity condition (22) up to order (kr) 2 then requires that that the 
fractional velocity perturbations for the scalar field and the radiation should be equal: v r ~ 
5<p/<j)'. Expressing the radiation velocity in terms of 5<p, the momentum constraint (last 
equation in (20)) then yields 

1 (^±m , (25) 




where p<f, = |0' 2 a~ 2 

The above equations may be used to determine the leading terms in an expansion in 
\kr\ of all the quantities of interest about the singularity. In order to compare with Ref. 
5, we shall choose to parameterize the expansions in terms of the parameters describing 



the comoving energy density perturbation, e m = —\Ti 2 k 2 Q, which has the following series 
expansion about r = 0: 



e m = e D(r) + e 2 E(r), (26) 
where eo and e 2 are arbitrary constants, and 

D(t) = 1 - 2r 4 r - ^Vm| kr\ + ..., 

E(t) = r 2 + . . . . (27) 
For adiabatic perturbations, we obtain 

5 * = 60 ("4^ " l lnlkTl + \ - I p) + £2 4^ + ° (T ' Tln|A;r|) 
V4, = e _ r 4 r )) + °( T i r\n\kr\), 

5r = ^ + 0(r 2 ,r 2 ln|A;r|), 
v r = v<p + 0(t, r\n\kr\), 

* = «» (~sh + S lll|fcT| + f I) " £2 sF + 0{T ' Tln| * T|) ' 

(S(j>) ( 3 „ „ , 1 , ,, , 1 13r?\ 1 , ,, ,, 



C 4 ,m = -^e 2 + e (^(A; 2 + 16r 2 ) + \ In |fcr|) + 0(r, r\n\kr\), ( 28) 



where C4,m is the curvature perturbation on comoving slices introduced by Mukhanov. 32 

In an expanding Universe the adiabatic growing mode corresponds to a curvature per- 
turbation, conveniently parameterized by C^m- The decaying mode perturbation is really a 
local time delay since the big bang, to which (^m is insensitive but $ is not. As detailed 
in Ref. 5, in a contracting Universe these modes switch roles so that the time delay mode 
is the growing perturbation and the curvature perturbation is the decaying perturbation as 
one approaches the big crunch. 

The perturbations generated in the ekpyrotic/cyclic scenarios consist of growing mode 
scale- invariant perturbations in the incoming state with no decaying mode component. 
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These perturbations are parametrized by eo/k 2 having a scale invariant spectrum, and since 
there is no decaying mode, (^m is zero on long wavelengths. After the collision, from the 
four-dimensional effective theory view the universe is expanding. Now, the growing mode 
perturbation is proportional to the long wavelength part of C,a,m- The key question is whether 
with our five-dimensional prescription matches the growing mode in the incoming state onto 
the growing mode in the outgoing state, parameterized by (^m, with nonzero amplitude. 
For this to occur, the long wavelength piece of (^m must jump across the bounce. We shall 
see below that this indeed occurs. 

IV. PROPAGATION OF GRAVITATIONAL PERTURBATIONS IN A 
COLLISION OF TENSIONLESS BRANES 

In this section, we consider the propagation of metric perturbations through a collision of 
tensionless branes where the background space-time is precisely M. /Z 2 x R 3 . The analysis 
follows closely Section II, which considered the propagation of generic scalar fields in this 
same background. The results here are essential to our analysis for the physically relevant 
case of colliding branes with tension (Section V) since our approach is based on finding 
a gauge where the propagation of metric perturbations through the bounce is as close as 
possible to the case for fixed tensionless branes. 

In this problem, it is simplest to choose coordinates in which the branes remain at fixed 
locations and all the fluctuations in the geometry are accounted for by the bulk metric 
perturbations. Recall that, ignoring gravity, the background metric is 

ds 2 = -dt 2 + t 2 dy 2 + dx 2 , (29) 

but with y identified under translations y — > y + 2y , and the reflection y — > 2y — y. The 
orbifold fixed points located at y — ±y /2 are the trajectories of two tensionless orbifold 
branes. In Section II we considered matching a scalar field across the singularity in this 
space-time 8 and now we generalize the methods considered there to the case of gravitational 
waves. 

A gravitational wave in five dimensions has five independent propagating components. 
If the y dependence may be ignored these five components split up in synchronous gauge 
into tensor (Sgij), vector (5g iy ) and scalar (Sg yy ) components, possessing two, two and one 
propagating degree of freedom respectively. As usual in four dimensional cosmological per- 
turbation theory the most interesting piece is the scalar as this transforms nontrivially under 
coordinate transformations and couples to the matter density perturbations. The tensor 
pieces are especially simple since they are trivially gauge invariant and decouple from the 
matter. Finally, the vector pieces only couple to the curl component of the matter velocities 
and not to the matter density perturbation. They require a separate analysis which will 
not be given here. Furthermore, in our setup the vector modes are naturally projected out 
because 5gi y must be odd under the Z 2 - Hence the vector modes must vanish on the branes, 
and this is why there are no vector degrees of freedom in the four-dimensional effective 
theory. 

We shall, therefore, need only to consider the scalar sector in what follows. The form we 
take for the five- dimensional cosmological background metric is 
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ds 2 = n 2 (t, y){-dt 2 + t 2 dy 2 ) + b 2 (t, y)5 ij dx i dx j , (30) 

and we write the most general scalar metric perturbation about this as 

ds 2 = n 2 (t, y)(-(l + 2<$>)dt 2 - 2Wdtdy + t 2 (l - 2Y)dy 2 
-2Wiadx l dt + 2t 2 V i (3dydx i ) 

+b 2 (t,y)((l - 2^)5ij - 2V i V j x)dx i dx j . (31) 

For perturbations on Ai c x R 3 it is straightforward to find a gauge in which the metric 
takes the form 

ds 2 = (1 + \k 2 X ){~dt 2 + t 2 dy 2 ) + ((1 - \k 2 x)ki + 2k i k jX )dx i dx j , (32) 

and x satisfies a massless scalar equation of motion on M c x R 3 . To be precise, the gauge 
is 

a = (3 = 0, r = $ - - A; 2 x, 

2 1 

$ = g^X, * = ^k 2 X , 

W = 0. (33) 
Notice that the non-zero variables can all be related to X according to 

(T,^) = (- 2 -,+ 2 -, + 1 -)k 2 x. (34) 

We shall, henceforth, refer to these as the 'Milne ratio conditions.' Furthermore, imposing 
the Z 2 symmetry, we obtain Neumann boundary conditions on X , 

x'(y±) = 0, (35) 

where y± = ±yo/2 are the locations of the two Z 2 fixed points. 

In the model space-time, the lowest energy mode for X is y-independent and has the 
asymptotic form 

X (t,y)=Q + P\n\kt\, (36) 

with Q and P being arbitrary constants, just like the case of scalar fields in Section II. Our 
matching proposal for all the perturbation modes is then simply the analogue of the scalar 
field rule given in Section II, namely 

Qout — —Qin + 2(7 — ln2)Pj„, Pout — Pin- (37) 

These relations are sufficient to determine the metric fluctuations after the bounce. 

In later applications, we are only interested in the long-wavelength part of the spectrum, 
and, for the cases of interest, P is suppressed by k 2 compared to Q. As a result, we obtain 
the approximate matching rule 

Qout Qin) Pout Pin- (38) 

The key conditions (33) through (35) are satisfied precisely for all time in a compactified 
Milne mod Z 2 background. When tension is added to the brane and the bulk is warped, 
our approach is to find a gauge which takes us as close as possible to these conditions in the 
limit as t tends to zero, where the same matching rule may then be applied. 
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V. 5D COSMOLOGICAL PERTURBATIONS FOR BRANES WITH TENSION IN 

A WARPED BACKGROUND 



Our strategy for computing propagation of perturbations when the branes are dynamical 
and have tension (so the bulk is warped) is conceptually simple: 

1. We use the four-dimensional effective (moduli) theory described in Section III to pro- 
vide boundary data for the five-dimensional bulk fields. In particular, we will be 
interested in the case where a nearly scale-invariant perturbations has been gener- 
ated well before the bounce when the four- dimensional effective theory is an excellent 
approximation, as occurs in ekpyrotic and cyclic models. 

2. In the five-dimensional theory, we find a gauge which approaches the Milne conditions 
(33) through (35) as t — > 0. In the gauge, the perturbation variables satisfy the 
massless scalar field equations of motion. 

3. We use the conditions in (38) to propagate all perturbation variables through the 
collision. 

4. We match onto the four-dimensional (moduli) theory to determine the cosmological 
results for long wavelength perturbations. 

One might worry that the four- dimensional effective theory we use to predict the bound- 
ary data for five- dimensional general relativity breaks down close to the bounce. However, 
there are reasons to expect the effective theory remains accurate as an approximation to 
general relativity even at small times. First, in Kaluza-Klein theory, the effective four- 
dimensional theory is a consistent truncation and hence provides exact solutions of the 
five-dimensional theory even in situations of strong curvature and anisotropy. In our case, 
as the branes come close, the warp factor should become irrelevant so that the Kaluza-Klein 
picture should become more and more valid. Second, in the approach to the singularity in 
general relativity 45 (based on the classic BKL work 46 ), the decomposition of fields according 
to dimensional reduction does correctly predict the asymptotics of the solutions in the limit 
as t — > 0. This suggests that the effective field theory indeed captures the correct behavior of 
full five- dimensional gravity near the singularity. In our detailed study of the linearized the- 
ory, we shall find a remarkable consistency between the predictions of the four-dimensional 
effective theory near t — and the full five- dimensional cosmological perturbation equations, 
and these consistency checks are the main justification for our use of the effective theory 
all the way to the brane collision. Of course, the use of five-dimensional general relativity 
near the singularity may itself be doubted since stringy corrections may be large there. But 
this objection can only be addressed in a detailed calculation within a string or M-theory 
context, which is beyond the scope of the present paper. 

We first infer the boundary geometry in longitudinal gauge (Section V.A) for which there 
is a simple and precise correspondence between the four- and five-dimensional perturbations 
and both are completely gauge fixed (see also Appendix 2). However, in this gauge the 
metric perturbations diverge much more rapidly (as 1/t 2 ) than a massless scalar near t — 0. 
We shall need to transform to a gauge where a) all the components of the metric are only 
logarithmically divergent and b) in which the components of the metric are in the same 
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ratios and obey the same boundary conditions asymptotically as t — > 0, as for the perturbed 
model spacetime with two tensionless branes in Ai° /Z 2 x -R 3 (Section V.C). In this gauge we 
can treat the components as massless fields and match across the singularity as in Section 
II (Section V.D). 

We wish to emphasize that the choice of gauge we are making is fully five-dimensional 
and is quite unlike that usually made in four-dimensional cosmology for several reasons. In 
four- dimensional cosmology, the matter present is often used to define a gauge - for example 
one may choose gauges in which the total density or velocity perturbation is zero. However 
in the five-dimensional bulk there is never any matter present, just the cosmo logical term 
which is constant and, therefore, does not define any preferred time-slicing. One might 
choose surfaces of constant extrinsic curvature, but these are not in any way preferred by 
the physics involved. Instead, our approach focuses on the asymptotic geometry near t — 0, 
and identifying it with the model space-time A4/Z 2 x R 3 . In addition to approximating 
the model space-time, it is essential that, for the same gauge choice, the brane collision be 
simultaneous at all x, so that the t = CT and t = + surfaces physically coincide. We shall 
show that our gauge choice satisfies this latter criterion, but the standard four- dimensional 
gauge choices, e.g., constant density or velocity gauges, do not. 



A. Longitudinal gauge moduli predictions 

In this section we wish to use the four-dimensional effective (moduli) theory discussed 
in Section III to infer the boundary data for the five dimensional bulk perturbations. In 
any four-dimensional gauge, the four- dimensional metric perturbation and scalar field 
perturbation 5<p determine the induced metric perturbations on the branes (in a related but 
not equivalent gauge) via the formulae (10): 

^ = V + 2 ( lnfi ±),<^0fiV> (39) 

where VL + = cosh(0/ \^6) and = — sinh(0/ v^6) and the metric perturbations are fractional 
i.e. 8g„ v = a 2 /^, 8g±„ = b\h± v . 

This formula is particularly easy to use in five- dimensional longitudinal gauge. (Our 
definition follows that of Ref. 50, where many useful formulae are given.) This gauge may 
always be chosen, and it is completely gauge fixed as we explain in Appendix 2. In this 
gauge the five- dimensional metric takes the form 

ds 2 = n 2 (t, y)(-(l + 2§ L )dt 2 - 2W L dtdy + t 2 (l - 2T L )dy 2 ) 
+b 2 (t,y)((l-2^ L )5 tJ )dx l dx^ (40) 

Furthermore, as explained in Appendix 2, in the absence of anisotropic stresses the brane 
trajectories are unperturbed in this gauge. An immediate consequence is that the four- 
dimensional longitudinal gauge scalar perturbation variables $± and \I/± describing pertur- 
bations of the induced geometry on each brane 

ds% = 6|(r±)(-(l + 2<& ± )dT 2 ± + (1 - 2V±)dx 2 ), (41) 

are precisely the boundary values of the five-dimensional longitudinal gauge perturbations 
$ ± = $l(j/±) and ty± = ^l(?/±)- Using (39) and (41), we find for the induced perturbations 
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$+ = $4 + —f= tanh(0/v / 6)50, 
V6 

^ + = $ 4 - ^tanh(0/v / 6)50, 

= $ 4 + -^= cath((j>/V6)S(j>, 
V6 

^_ = $ 4 - -= coth(0/ v / 6)50. 
V6 



(42) 



One subtlety in utilizing these formulas is that if $ 4 and 5(f) are expressed as functions of 
four-dimensional conformal time, then they give the correct predictions for $± and \I/± on 
the branes in terms of the conformal time t± on each brane. However, when we use them 
as boundary values of the five- dimensional metric it will be necessary to consider all the 
perturbation variables as functions of the five- dimensional time t entering in the background 
metric (30). The brane conformal times may be expressed in terms of t by integrating, 



Jo oft, ?/+) 



where q = b/n. So for example the boundary value of the bulk metric perturbation on 
the positive tension brane is given explicitly by 

$ L (t, y+ ) = $ 4 (| q{t,y + )' l dt) 

+ tanh(0( J q(t, y+^dt) / 'V6)S<j>( J q(t, y+^dt), (44) 

where y + is the location of the positive tension brane. As noted, in this gauge even when 
we include perturbations the branes are static and the Israel matching conditions are easily 
found to be 

b' L 

t(v±) = ±-^ nt P±, 

b 6 

%±) = ±^nt(p± + p±), (45) 
q z 

for the background solution and 

K(y±) = \w L =F \nt{8pt - r LP± ), 

= + g- t W L =F -^nt(5pt - T LP± ) =F -nt(5pt - T L p ± ), 

W L (y ± ) = ±—( P± + p ± )vt, (46) 
n 

for the perturbations, where the right hand sides are all evaluated at y — y±, the locations 
of the positive and negative tension branes. (From now on prime shall denote d/dy and dot 
shall denote d/dt.) We can re-express Wl on the branes as 
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W L (y ± ) = (q 2 )'vt, 



(47) 



where v L is the longitudinal gauge velocity perturbation of the matter on each brane. From 
this one sees for example that for empty branes, Wl vanishes on the branes. 

As long as the bulk matter is isotropic, as it is in our case, the Einstein equations lead 
to a constraint which may be written 

G\ - G 2 2 = 0. (48) 

In longitudinal gauge this reads 50 

r L = $ L - * L , (49) 

everywhere in the bulk. This is the five-dimensional analogue of the well known four- 
dimensional no-shear condition $ = ^ in longitudinal gauge. Equation (49) serves to define 
on the branes in longitudinal gauge. Consequently we have sufficient boundary data 
for all the components of the five- dimensional metric in this gauge. We can then perform 
an arbitrary five-dimensional diffeomorphism to infer the boundary data in any gauge we 
choose. Equivalently, equation (49) may be interpreted as a condition in any gauge by using 
the gauge invariant variables defined in Appendix 2. 



B. Stress energy conservation 

In this paper, we consider perturbations in the 'in' state which may be described as 
local fluctuations in a single scalar field <fi representing the inter-brane separation. We 
are interested in long wavelength modes which are completely frozen-in during the collision 
event. Hence the local processes describing the production of radiation at the bounce should 
be identical at each x, and in the usual sense employed in cosmology, described in Section 
HID, the perturbations should be 'adiabatic'. 

As is well known, the conservation of stress energy leads to powerful constraints on 
adiabatic density perturbations, in particular implying that the amplitude of the growing 
mode perturbation cannot be altered on super-horizon scales. In this section we discuss this 
constraint and show how it implies the spatial curvature of comoving (or constant energy 
density) slices is conserved on large scales both for the brane geometries and for the four- 
dimensional effective theory. We shall restrict ourselves to considering only radiation on each 
brane. This considerably simplifies the analysis because when the matter on each brane is 
conformally invariant, as explained above, in the four-dimensional effective theory the scalar 
field decouples from the matter and can be treated as an independent fluid. 

First we need to generalize the usual notion of adiabaticity to deal with perturbations 
in the radiation densities on each brane. As mentioned above, radiation couples to the scale 
factor Q + (<f))a on the positive tension brane and f2_(0)a on the negative tension brane with 
notation as in the previous section. Conservation of the radiation density on each brane 
reduces at long wavelengths to dp± = — 4p±d\na — 4p±d\nQ,±. Likewise we have for the 
radiation density in the four dimensional effective theory dp A = — 4p 4 rflna. Hence solving 
for 5 In a as in Section HID, we infer the adiabaticity condition for radiation on the branes 
to be 
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6± = 5 A -4(lnfi±)^50. 



(50) 



The equation for conservation of energy in four dimensions can be written in the form 50 

Cb = \k 2 v L , (51) 

where 

Cfl = *-3(TW' <52) 

is the gauge-invariant variable measuring the spatial curvature perturbation on constant 
density hypersurfaces, as originally defined by Bardeen. 49 The quantity vl is the gauge 
invariant scalar velocity potential, equal to the velocity potential in longitudinal gauge (so 
that VjfL is the scalar part of the velocity perturbation). 

At long wavelengths k — > 0, equation (51) implies that (b is conserved, provided the 
velocity perturbation does not grow with scale. This property is very powerful since it 
means that under most circumstances, as long as modes remain outside the horizon (b can 
be trivially extrapolated from the early to the late Universe, where it gives the amplitude 
of the growing mode adiabatic density perturbation, the main quantity of observational 
interest today. 

The above definition (52) applies equally on each brane and in the four-dimensional 
effective theory, provided the terms on the right hand side are appropriately interpreted. 
On the branes, we have 

Cb,± = *± - \S ±: (53) 

where 5± are the fractional perturbations in the radiation densities on each brane, and \I/± 
is the perturbation in the brane spatial metric. Using (42), written as 

# ± = ^ 4 -(lnft ± )>/>, (54) 

and the adiabaticity condition (50) we see that the four-dimensional effective value of 
Bardeen's variable, Cb,4 = $4 — \8± is in fact identical to Cb,± on long wavelengths. 

Our final result will in fact more naturally emerge in terms of another gauge invari- 
ant variable, the curvature perturbation on comoving slices, emphasized by Mukhanov and 
others. 32 This is defined as 



{ M = * + Hv, (55) 

with v the velocity potential and 7i = d\na(r) j 'dr the conformal Hubble constant. Again 
this may be interpreted on either brane or in the four-dimensional effective theory. But 
adiabaticity requires that the fluid velocities be identical on long wavelengths for each fluid 
component. Therefore we must have v± = V4 = 5<f)/<f)' (from (23). This is also seen to be 
consistent with (51) and the equality of the Bardeen variables (b,± = (b,4 which we have 
just shown. 
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The scale factors on each brane are related to the four-dimensional effective scale factor 
via b± = Q±a. Recalling that the conformal times on the branes are the same as that in the 
effective theory, we have 7i 4 = TC± — (In ^±),</>0, T - Using v 4 = v ± = = 6<f>/(<j> iT ) we find 

Cm a = ^4 + n 4 v 4 = ^ 4 + n±v± - (lnOtX^ty 

= * ± + H±v^ (56) 

which is just Cm,±- So for adiabatic perturbations and at long wavelengths, the comoving 
curvature perturbations on the branes are both equal to that in the four- dimensional effective 
theory. As is well known, the latter is conserved for for adiabatic perturbations at long 
wavelengths. It follows that away from the bounce, Cm,± are both conserved as well. As we 
discussed in the introduction, and will detail below, this does not imply they are conserved 
across the bounce. 

We will use (56) below, but we should point out one minor subtlety. We shall be perform- 
ing all our calculations in five-dimensional time t, not four-dimensional conformal time. The 
velocity is not a scalar under coordinate transformations, and we shall need to multiply 

by a factor of q when we re-interpret equation (56) in terms of the five-dimensional time 

t. 



C. Transformation to Milne gauge 

Our philosophy is to evolve cosmological perturbations through the bounce in a 'Milne 
gauge' where they behave as closely as possible to gravitational waves on Ai c /Z 2 x R 3 , as 
described in Section IV. Then, we can use the same matching conditions (38) to determine 
the perturbation spectrum after the bounce. 

The Milne gauge we use is chosen to match the gauge choice (33) in Section IV up to 
corrections of order t and tin \kt\ due to the finite brane tension, radiation densities and the 
warp factor. We still have enough coordinate freedom to set three linear combinations of 
the metric perturbations equal to zero for all t, and we choose 

a = (3 = 0, r = $ - - k 2 X - (57) 

A remarkable feature of this choice is that the constraint equation (49) implies that X obeys 
the equation for a massless scalar field on the unperturbed background for all times. From 
(57) and (99) in Appendix 2, we find 

vy2 Id d X ,ld 2 X bdx^3b'd X k 2 b 2 

v x = -tdt {t m ] + + ¥ ~b T y -^ x = - (58) 

This result is remarkable in that it is independent of the precise details of the back- 
ground bulk geometry and the form of the stress energy in the bulk, assuming only that no 
anisotropic stresses are present. 

The remaining gauge freedom is of the form x M — > x M + £ M where 

b 2 ■ b 2 

e = V> e = (59) 
n n t 
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provided that £ s also satisfies a massless scalar field equation 



v 2 e = o. 



(60) 



Since x transforms as \ — > x + £ s > an d X is zero i n longitudinal gauge, it follows that x i n the 
gauge we use is, in fact, precisely value of the spatial coordinate transformation £ s needed 
to get to a gauge satisfying (57) from five- dimensional longitudinal gauge. Furthermore, £' 
and £ y may be inferred from x — £, s via (59). 

To completely fix the gauge within the family specified by (57), we need to specify 
boundary conditions for the field x on the two branes, and initial conditions on some space- 
like surface. As a first guess, one might consider choosing to fix the gauge by specifying 
Neumann boundary conditions on the branes (i.e. ^{t, y±) = 0) for all time, as in Section IV. 
One can easily prove that in this gauge, as in longitudinal gauge, the brane trajectories are 
unperturbed. This follows from the formula (59) upon setting f = x as noted above. This 
is very important: it follows that in this Neumann gauge the brane collision is simultaneous 
and occurs at precisely t — for all x. Furthermore the Neumann gauge x?(t,y±) = for 
all t is a good gauge in the sense that none of the metric components diverge worse than 
logarithmically. 

However, it turns out that setting x?(t,y±) = f° r a U time is too strong a condition. 
One cannot choose Neumann gauge for all time and also have 



consistent with the behavior in the model space-time (33) at leading order in t and tin \kt\. 
The resolution is simple: we need to perform a small gauge transformation away from 
Neumann gauge in which we maintain only the asymptotic vanishing of the proper normal 
derivative of x as t tends to zero, i.e. we impose that 



on the two branes. With this choice we are able to impose all of the conditions in (63) as well 
as (57). This small gauge transformation away from Neumann gauge shifts the locations of 
the branes, y±, but only by a finite amount. As discussed in Appendix 4, this means that 
the rapidities of the branes are perturbed in our chosen gauge, but the collision event is still 
simultaneous. 

Our reason for expecting that we can choose a gauge specified by (63) and (64) is that 
when the branes approach the warp factor should become increasingly irrelevant and the real 
background space-time should asymptotically approach the model space-time M./Z2 x R 3 . 
We expect the low energy modes we are interested in to behave as the lowest Kaluza-Klein 
modes in this limit, i.e. becoming independent of y. Within the class of gauges specified 
by (57), we shall indeed see that there are solutions for the perturbations in which all 
the perturbation components behave like Q + Pln|/c£| as t tends to zero. The Milne ratio 
condition (34) turns out to be automatically satisfied by the coefficients of the logarithms. 



W = + 0(t, t\n\kt\) 
$ = h 2 x + 0(t,t\n\kt\) 

^= l -k 2 X + 0{t ) t\n\kt\), 



(61) 
(62) 



(63) 



n-VY(y ± ) = O + 0(MlnN), 



(64) 
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Fixing the constant terms to be in the Milne ratios further fixes the gauge up to a residual 
two-parameter family and imposing asympotically Neumann boundary conditions (64) on 
both branes then completely fixes the gauge. 

Imposing asymptotically Neumann boundary conditions turns out to have various other 
natural consequences. For example in this gauge, all the metric perturbation components 
possess identical asymptotic behavior (i.e. constant and logarithmic terms) on the two 
branes, as t tends to zero, consistent with their behavior as a lowest Kaluza-Klein mode. 
Furthermore, there is a simple geometrical consequence of this choice which we explain 
in Appendix 4, namely that the in this gauge the perturbations to the embedding (T, Y) 
coordinates of the brane collision event actually vanish so the branes collide at precisely the 
background values of T and Y. 

The non-zero perturbations in our chosen class of gauges are $, W and x along with 
T which is fixed by the gauge choice (57). All the gauge freedom is contained in the solution 
for x- To see this we note that if we know the solutions for x we can immediately infer <3>, \l> 
and W from the values in longitudinal gauge via the formulae from Appendix 2, (99), which 
with (57) imply 

/ 6 -\ ™ 2 • n> f <fx\ 

$ = $ L - _ _ g 2 x + _ *A 

n n t z 



W = W L -(q 2 x)'-t 2 (^f 



Here as above, q = b/n. 

Our goal then is simply to determine x to sufficient order in t to be able to compute 
all the other components from (65). As we have already explained, in our chosen class 
of gauges x satisfies the massless scalar equation (58) at all times. To specify a complete 
solution we need to specify both Cauchy data on some constant t hypersurface between the 
two brane worldsheets, plus boundary conditions on the two branes. The boundary data 
will be obtained from the four-dimensional effective theory, and we make the conjecture 
that the bulk solution which is consistent with these data will behave near t = like a 
Kaluza-Klein zero mode on M c / Z 2 x R 3 , which is to say that the perturbations should be 
independent of y as t tends to zero. In practice this means we will look for a solution which 
is asymptotically of the form x — Q + -Phi \ kt\, independent of y. This assumption formally 
provides the Cauchy data once we determine Q and P (see below). 

At higher orders in t, we shall allow for arbitrary Neumann boundary conditions, which 
we shall parameterize as 

X\V±) = \4t 2 + 0{t\tHn\kt\). (66) 

As explained above, we shall adjust the coefficients af to obtain the correct Milne ratios. 
Note that there can be no 0(1) term since we are assuming that x is asymptotically of the 
form x — Q + -Pint, independent of y, and the 0(t) term is prohibited by our condition 
(64). In principle we could also include tint and t 2 lnt terms but we shall find that the 
Milne ratio conditions (34) are sufficient to rule these terms out. 
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The form of the series expansion for x, implied by its equation of motion (58), is 



X (t, y) = (Q + (fi(y) + ci coshy + c 2 sinhy)* + f 2 (y)t 2 /2 + 0(t 3 )) 
+P\n\kt\(l-^kH 2 + 0(t 3 )), (67) 

where /i(y) and /2(y) are two functions of y that are obtained as solutions of second order 
differential equations in y with boundary conditions derived from (66). We choose to define 
fi so that f[(y±) = 0. Therefore if x satisfies the asymptotically Neumann condition (64) 
on both branes, we must have c\ = c 2 = 0. A geometrical interpretation of this condition is 
explained in Appendix 4. 

Using the expressions for b(t,y) given in Appendix 3, equation (123), in the equation of 
motion (58), for x, we find at order the following differential equation must be satisfied 

by/i(y): 

fi ~h- 2L ^ nhyQ ((6 + r + L 2 ) cosh(y + |) - (6 - r_L 2 ) cosh(y - |)) = 0, (68) 

A similar equation for f 2 is found at order t°. The solutions are messy in general but simpler 
when no radiation is present, for example in this case we have 

fi(y) = — — — i~r fvcoshy — (1 + — tanh— ) sinh-u^) . (69) 
J yy > 2Lcosh(?/ /2) V V 2 2 ' y J y J 

By substituting (67) into (58) and imposing the boundary conditions (66) at each order, 
the solution for \ U P to t 3 corrections is completely determined in terms of the four constants 
in total: Q, P, and af. From this solution for x, equations (65) then determine all the other 
components of the metric perturbations at leading order in t, on each brane. 

Let us start by determining the spatial curvature perturbation \l/ on each brane. From 
(56) and (65) we find 

b b' a V 

^ = a,M + Q l (QX-v (f) )--^. (70) 

We require that ^ be only logarithmically divergent. Since from (28) we have that v$ = 
3e /(4A; 2 r) + 0(1), diverging as t^ 1 as t — > 0, we see from (67) and the expressions for the 
background metric functions in Appendix 3 that only \ can cancel that divergence, which 
requires that 

p = w (71) 

This condition ensures that the curvature VP in our gauge and the comoving curvature ^4,m 
in the four dimensional effective theory only differ by a constant at leading order in t. 
However, it shall be very important that the constant is nonzero. As we shall see, the 
constant represents the time delay between the two time-slicings, and it is the key to why 
jumps across the singularity. 
We shall now show that it is possible to choose the three remaining gauge constants Q 
and af so that the metric takes the canonical Milne gauge form asymptotically as t tends to 
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zero. First, in this gauge all the metric perturbations behave as Q + Pln|/ct|, as t tends to 
zero, but with different constants Q and P for each component. Substituting (67) and (71) 
into (65), with ^ ^ given from (56), $£, given from (44) and Wl given from (47), one finds 
that the logarithmic terms are actually all in the correct Milne ratios (34), and also that 
W vanishes to leading order, independently of the undetermined constants. Furthermore, 
the logarithmic terms obey 3>(y+) — $(y-) = and ^(y+) — ^(y-) = 0, consistent with our 
assumption that the Kaluza Klein zero mode dominates. 

The gauge constants Q, af and af do, however, affect the t-independent constant terms 
in each metric perturbation component. Two of the constants are fixed once one sets the 
constant terms in $ and ^ to their Milne ratio values (2/3)k 2 x and (l/3)A; 2 x- We also 
want to ensure that all components of the metric perturbations behave asymptotically like 
a Kaluza-Klein zero-mode, becoming independent of y as t tends to zero. We check this by 
comparing the values of <E>, \& and W on the two branes. The difference of ^ on the two 
branes turns out to be independent of the choice of the gauge constants as t tends to zero, 

=0(p 2 ± L 2 ). (72) 

Since the moduli space approximation was derived neglecting p\L 4 corrections, to the order 
we can trust the calculation, \1> is equal at the two brane locations. The difference of $ 
on the two branes is not automatically zero at leading order however, and setting it zero 
provides the additional equation needed to determine the third constant. 

The result of these calculations is that the solution for x U P to 0(t 3 ) and the leading 
order behavior of the other components of the metric is completely determined. Explicitly 
we find 47 

m = (4Mt) + 32fc°LW(y /2) " Sin%0)) " ^ (r+ " r - )( " % ° + Sin%0) ) 

+0(p 2 ± L 2 ,t,t\n\kt\). (73) 

Since ^ is one of the variables which we match in our chosen gauge, it follows that our 
prescription is quite different to matching the comoving curvature perturbation (^m four- 
dimensional effective theory. As we shall explain, the additional terms in (73) allow the 
propagation of growing mode perturbations across the singularity. 



D. Matching Proposal 

The requirement that around the collision event the geometry looks locally like Ai /Z 2 x 
R 3 has completely fixed the gauge in the incoming and outgoing states. As elaborated in 
Appendix 4, the asymptotically Neumann boundary condition (64) further ensures that the 
collision event is simultaneous in our gauge, an essential property for matching perturbations 
since the space-like surfaces defined by t — > + and t — > 0~ then physically coincide. 

Furthermore as we have discussed this gauge is special in that the induced geometry on 
each brane is asymptotically the same at collision. In general if a brane is moving, the values 
of the bulk perturbations $ and * evaluated on the branes differ from the induced values 
$ ± and The differences are given by 
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WJ_/ n 2 t 2 n A 
b 2 b' 

* ± - = +^ T x: 

Since n'/n oc t and oc t as t — > 0, in the presence of matter on the branes, if we 
make the requirement that the metrics on each brane are asymptotically identical this fixes 
X' = + 0(t 2 ), which is what we have required. Physically this seems a natural choice of 
gauge because when two ordinary branes collide, the induced geometries are identical at the 
collision moment. This interpretation is also consistent with the predictions from the four- 
dimensional effective theory where the brane metrics are give by g+ = (cosh(0/v^6)) 2 5' /tl / 
and g~ v = (— smh((j)/^/6)) 2 g^ v . Since the collision corresponds to — > — oo, in the limit 
there is no difference between the two conformal factors, and the brane geometries appear 
identical. 

So in the gauge we have fixed by the requirement that the perturbations behave asymp- 
totically like those on the model space-time Ai/Z 2 x R 3 , the collision event is synchronous, 
the Milne ratio conditions are satisfied, the boundary conditions are asymptotically Neu- 
mann, and the geometries on each brane are asymptotically the same both before and after 
the collision. Our matching proposal amounts to relating the geometry on these chosen time 
slices across the collision. We believe that these are sufficiently desirable properties to justify 
this as the natural gauge in which to perform the matching, and from now on we shall take 
this to be our complete gauge fixed matching gauge. 

Let us now return to our final formula (73) to infer its meaning in the context of ekpy- 
rotic and cyclic models. In those scenarios 5 the quantity eo/k 2 has an approximately scale- 
invariant long wavelength spectrum in the incoming state. The first point to make is that 
even in an in-state with no radiation present the dimensionless curvature perturbation on 
spatial slices \I/ in our gauge has a scale invariant spectrum, since 

Recall that y is the relative rapidity and V in = tanh(y ) is the incoming relative velocity 
between the two branes. Then, at small velocities this gives 

* = ~ l±¥b v ™- (75) 

We may interpret this geometrically as follows. In the absence of radiation there is no 
real meaning to the curvature perturbation on the branes but if we imagine that there is a 
small density of radiation coming in, and the perturbations are adiabatic, we can infer the 
comoving curvature perturbation on the brane, C±,m, so (75) becomes for long wavelengths 

(76) 

Since C±,m and \I/ are the spatial curvature perturbations of the branes as respectively 
measured in the comoving timeslicing and in our chosen timeslicing (in which the collision is 
at t — 0), it must be that the additional piece arises from a time translation between the two 
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gauges. That this is so is verified when one traces back the origin of this term to the second 
term in equation (70). As explained in the introduction, comoving gauge (or equivalent 
constant energy density gauge) are bad gauges to match in because the brane collision is not 
simultaneous in those gauges. Since our prescription is to propagate \I/ across the collision, 
the jump in (± is due to the time delay occurring between the collision-synchronous surfaces 
in our gauge, and those of the comoving/ const ant density surfaces. The key to our result is 
that in the comoving or constant density gauge the time delay between t = and the actual 
brane collision event has a scale invariant spectrum. 

In fact, using (28) and (73) to find Q and P before and after the bounce for all compo- 
nents of the metric perturbations and matching according to the rule given in equation (38) 
results in (^m inheriting two separate scale-invariant long wavelength contributions in the 
post-singularity state. The first occurs as a direct consequence of the sign change in (38), 
and is independent of the amount of radiation generated at the singularity. The second is 
proportional to the difference in the densities of the radiation on the two branes. At leading 
order in velocities we have 

AC4,M = l^VL + O " W ^ + 0(r±V "' V5L " ^ (77) 

where V in and V out are the relative velocities of the branes before and after collision. Note 
that since P oc eo, matching P is in fact equivalent to matching eo across the collision 
as proposed in Ref. 5. In terms of four-dimensional parameters defined in Section III.C 
including given in (18) defining the abundance of the radiation and the fractional density 
mismatch / defined in (17), we find again at leading order in velocities 

a a _ 3 e o (v a , t/4 \ 3e fr A V* ut 

AU > M ~ umJ [Vm + v ° ut) ~ i6F 1 ■ (78) 

This is our final result, relevant to tracking perturbations across the singularity in the 
ekpyrotic and cyclic models. We see it consists of two essentially independent terms. The 
first is proportional the radiation density mismatch on the two branes after collision. Note 
that just such a mismatch (with more radiation on the negative tension brane) was required 
in order to enable the cyclic solution of Ref. 2 to work. The second term exists however even 
in the limit of no radiation generated on the branes. As we have noted above, it is nonzero 
even if Vi n = V ou t, and it originates in the sign change of the parameter Q in our matching 
rule, which yields an arrow of time across the collision as explained in Section II. Going back 
to the original formula (73) in which we have not made the small velocity approximation, 
we note that both the radiation-dependent and radiation-independent terms possess a well 
defined limit as the brane collision becomes relativistic (as the rapidity yo — > oo), 



AC4,„«S ^^^^ • (79) 




Recall, we need the radiation densities on the branes to be much smaller than their tension 
i.e. r±L 2 < 1, in order that the four-dimensional effective theory be valid (Section III). 
Therefore in the high velocity limit, the radiation-independent term dominates. Conversely, 
from (78), in the low velocity limit (with (r+ — r )L 2 fixed) the radiation-dependent term 
dominates. 
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We should stress once more that the dependence upon parameters in (78) indicates its 
thoroughly five-dimensional origin. It cannot be expressed in purely four-dimensional terms. 
In previous work 5 with Khoury and Ovrut, two of us employed a more naive matching pre- 
scription framed entirely in terms of the four-dimensional effective theory. This prescription 
was based upon using the comoving energy density perturbation e m , which is finite at the 
singularity, as the matching variable. Unfortunately, since the differential equation govern- 
ing e m is singular at t — 0, the first time derivative of e m is not an independent quantity at 
the collision and hence could not be independently matched. Instead we proposed matching 
the second time derivative. This has the virtue of at least yielding a dimensionally correct 
result, but it is ambiguous since there are other choices of finite variables. Now we under- 
stand the source of the ambiguity better. There is simply not enough information present 
in the four-dimensional theory to fix the gauge. For that, the five-dimensional picture is 
essential as we have seen here. 

In summary, we have found that a spectrum of scale invariant, growing, long wavelength 
perturbations generally propagate across the singularity even in the limit when no radiation 
is produced. The radiation-independent contribution rests upon the sign change of Q in the 
matching rule (5). If radiation is produced at the bounce, then, for the long wavelength 
modes we are interested in, we believe it is reasonable to model the production of radiation 
as occurring suddenly, taking into account the conservation of energy and momentum as 
was done in Ref. 2. In this case, we find an additional contribution to the long wavelength 
scale-invariant perturbations emerging from the singularity, which is proportional to the 
difference in the radiation densities on the two branes. 

VI. CONCLUSIONS 

In this paper we have developed an unambiguous and, we believe, compelling rule for 
matching perturbations across the types of singularity encountered in the ekpyrotic and 
cyclic Universe scenarios. In the simplest realization of these scenarios, involving the collision 
of two Z 2 branes in a bulk with a negative cosmological constant, we have shown that 
the proposed rule leads unambiguously to a spectrum of scale invariant growing density 
perturbations in the ensuing hot big bang phase, even in the limit when only a small amount 
of radiation is produced at the collision. The result provides support for a key assumption 
of the ekpyrotic and cyclic models. 

We have dealt here only with the linear theory, treating the perturbations as free massless 
fields which we match across the singularity. This treatment clearly is not fully consistent 
since the perturbations are divergent at the singularity and nonlinear effects must become 
important there. However, there are reasons to expect that in the nonlinear theory, a similar 
matching rule will apply. In the linear theory, we have seen that the metric components 
typically behave as 1 + e\n\t\ as t — > 0. This is just the small e expansion of \t\ e , the generic 
behavior expected in the full nonlinear Kasner solutions of general relativity that describe 
the generic approach to a space-like singularity. The natural extension of our proposal to 
the nonlinear theory is, therefore, that we should match the Kasner exponents across the 
singularity. As in the linear theory, the canonical momenta associated with the three-metric 
are finite and our proposal amounts to matching them with a sign flip. But our matching 
proposal in (5) also reverses the long wavelength component of the constant term in the 
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metric perturbation. Generalizing to the nonlinear case, we may anticipate that when the 
metric tends to the Kasner form with spacelike components ~ e^|t| p , with Q and P of order 
e, these components will match to e~ Q |t| p in the outgoing state. If e is small as expected in 
the ekpyrotic/cyclic scenarios, nonlinear corrections will be of order e 2 and hence negligible. 

Strongly rupporting the idea of a local matching rule is the classic conjecture that in 
general relativity the behavior of the metric and fields becomes ultralocal in the approach 
to this type of singularity. 46 That is, the spatial derivatives become unimportant and the 
geometry at each point in space follows a homogeneous Friedmann-Robertson Walker evo- 
lution that just depends on local conditions. One might worry that contraction also leads to 
chaotic mixmaster behavior in which the universe moves from one kind of Kasner contrac- 
tion to another and the Kasner exponents change unpredictably. However, the existence of 
mixmaster behavior depends on the number and types of fields. We discuss elsewhere 48 how 
the mixmaster behavior can be naturally avoided in ekpyrotic and cyclic models. 

Finally, with a precise matching rule for propagating perturbations through the singu- 
larity in place, we believe that the cyclic and ekpyrotic models are now on firmer footing. 
A detailed study applying the above results to these cosmological scenarios will be given 
elsewhere. 52 
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APPENDIX 1: PROJECTED EINSTEIN EQUATION 

In Section III we derived the four- dimensional effective action for solutions with cosmo- 
logical symmetry and then used general co-ordinate invariance to infer the covariant moduli 
space action. While we have shown that this approach recovers the cosmological solutions 
perfectly at low densities, we have obtained the low energy effective action describing the 
general (asymmetrical) case by simply assuming locality and imposing covariance. While 
this is plausible it is important to check it explicitly. This has in fact been done in Refs. 
41,42,43 which further clarify the conditions under which the moduli approximation is valid. 
We shall compare the results of these works with those of the moduli space approach. 

We shall first show that the effective theory we have derived satisfies one non-trivial 
check. One way of formulating a low energy theory for the brane geometries is the so called 
Gauss-Codazzi formalism developed in Ref. 44. Here we take the five-dimensional Einstein 
equations and project them onto the brane to infer an equation for the brane geometry. One 
finds 
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where T^ u is the stress-energy on the brane, not including the tension. This looks like the 
four-dimensional Einstein equations except for two additional source terms. One contains 
stress energy squared terms, 

S»„ = \-TT, u - -T m T: + 1^(3T Q/3 T^ - T 2 ), (81) 

where T = T A A , and the second E^ v is obtained from projecting the 'electric' part of the 
bulk Weyl tensor onto the brane 

dx^ dx^ jj 
E^v = ~q^~q~^Eab } Eab = C acbdti n , (82) 

where n A is the normal to the brane. Note that by definition E^ v is symmetric. Since this 
term contains information about the second i y > derivatives of the bulk geometry we cannot 
calculate it in any purely four-dimensional way and so although the above equations strongly 
resemble Einstein's equations they are purely formal. However we can construct one purely 
four dimensional equation because E^ u does satisfy the exact condition 

E±n = 0. (83) 

The moduli space approximation only works in the limit in which the stress-energy of the 
matter on the brane is much smaller than the brane tension. This amounts to neglecting 
the T 2 terms in the above action leaving 



From now on we shall for convenience use units where M4 = {&nG)~ l l 2 is unity. As a 
consequence of the Bianchi identities it follows that in this 'low energy' approximation the 
following condition must be true. 

= 0. (85) 

Since E^ u is conserved and traceless it means that the influence of the bulk on the brane 
geometry is identical in form to that of the stress energy of a conformal field theory. If we 
look for a cosmological solution, the vanishing trace condition tells us that the only non-zero 
components of E£ are, E® = f(b) and E l - = — \f{b)5) where f(b) is an arbitrary function 
of the scale factor on the brane. In addition the condition that V nE± v = tells us that 
f(b) = C/b 4 and so the effect of this term is gravitationally indistinguishable from radiation, 
and it may be thought of as a dark radiation term. This is the import of Birkhoff 's theorem 
in the bulk, viewed from the brane. 

The moduli space approximation as we have developed it provides a precise prediction 
for E^ v . A non-trivial check on this approximation is that the predicted value of E^ v is 
traceless. This condition of tracelessness is built in at the start in the other formalisms, 41-43 
but is a nontrivial check of our approach. We can compute the trace by simply conformally 
transforming the trace of the Einstein equation in the four-dimensional effective theory. 
Writing the brane metrics as g^ u dx^ 'dx v = Q\g^ v dx^dx u we find 
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E% = -G^±T ± = R ± ±T ± 

= n ± 2 (R-^-v 2 n ± )±r ± 

= ft ± 2 (-T 4 + (V0) 2 (1 - 6 %i^) - ^%1^V 2 0) ± T± 

= ^ 2 (-T 4 -^^V 2 0)±T ± , (86) 

where T± = T±n and in the last step we have used Q + = cosh(0/y / 6), fi_ = — sinh(0/y / 6). 
Finally, making use of the equation of motion for the scalar field 

V 2 <P = -\(ni),,T + - 1 -({li), (P T-, (87) 

we find that 

= 0. (88) 

It is interesting to note that the intermediate steps in this calculation require that the 
conformal factors on the positive and negative tension branes are of the forms described 
above involving cosh(0/v^6) or sinh(0/v / 6)- 

In order to compute the projected Weyl curvature in general it is helpful to work at the 
level of the action. We start with the action for the four-dimensional effective theory 

S = J d 4 x^(^(R - (V0) 2 ) + Q 4 + C+ + (89) 

To get the action for the metric on the positive tension brane we simply perform the con- 
formal transformation, taking us out of Einstein frame 

s = J d 4 x^^n- 4 (^-(R + - 6^+v 2 ^; 1 - (v+0) 2 ) + n 4 + c+ + niC-), (90) 

then defining ^ = Q+ 2 and performing an integration by parts we obtain the following action 
for the metric on the positive tension brane 

S + = J d 4 xV^(\^R + - ^-^y (V+*) 2 ) + JC + + (1 - *) 2 £_), (91) 

and a similar calculation on the negative tension brane defining $ = Qz 2 gives 

S- = J d 4 xV^(\(*R- + (V_$) 2 ) + + (1 + $) 2 £_). (92) 

These results are in perfect agreement with the low energy approximation developed in 
Refs. 41,42 using a metric based approach and in Ref. 43 using the covariant curvature 
formalism. After deriving the equations of motion by varying these actions we can simply 
read off the predictions for the projected Weyl tensor on the positive tension brane as 
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1 (l-^) 2 

~ L(i-^ (v ^ v ^ - 1 6 ^^' w 

and on the negative tension brane 

1 fl + $) 2 
E\ = -TH V {1 + -)- [ ~^T% 

-^-(v^v;$-^v 2 $) 

+ l $(i + $) (v - $v:$ ~ \ 5 ^ v -^ ■ ^ 

A specially interesting limit of these equations is obtained by — > implying ^ — > 1 and 
$ — > oo which corresponds to the distance between the branes becoming infinite. Providing 
we can neglect the derivative terms we see that in this limit, matter on the positive tension 
brane couples to the brane geometry by means of the conventional four-dimensional Einstein 
equations, whereas the geometry on the negative tension brane is dominated by its coupling 
to matter on the positive tension brane, and will only start to look like conventional Einstein 
gravity if a 'stabilization' mechanism exists which freezes to a constant value. In the latter 
case, stress energy on each brane acts like a dark matter source for gravity on the other brane. 

These equations (94) describe matter interacting in an unconventional way with gravity, 
and yield a more complicated perturbation theory than usual. Our approach makes it 
clear that it is simpler to work with the effective four-dimensional theory in Einstein frame 



+ 



with a scalar field with a canonical kinetic term, and then simply to use the map g 
(cosh (0/ a/6)) 2 <? M „ and g~ v = (— sinh((j) / \^6)) 2 g ^ u to infer the brane geometries. The only 
sense in which this theory differs from conventional four-dimensional physics is that the 
different forms of matter couple non-minimally to gravity though the scalar field. 



APPENDIX 2: GAUGE INVARIANT VARIABLES 

As in four-dimensions the cosmological symmetry of the background metric allows us 
to find a set of gauge invariant variables, which facilitates the comparison of two differ- 
ent gauges. What the natural gauge invariant variables are depends on the form of the 
background and our definition closely follows but are not identical to those in Ref. 50. 

We begin with the background metric written in the form 

ds 2 = n 2 (t, y)(-dt 2 + t 2 dy 2 ) + b 2 (t, y)5 ij dx i dx j . (95) 

We shall only consider spatially flat cosmologies for simplicity but the generalization to 
closed and open universes is easy. The most general scalar metric perturbation can be 
written as 
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ds 2 = n 2 (-(l + 2<&)dt 2 - 2Wdtdy + t 2 (l - 2T)dy 2 ) 
-VS/iadxidt + 2t 2 V i f3dydx i 

+b 2 ((l - 2^)5ij - 2V i V j x)dx i dx j , (96) 

writing the perturbed metric as qab + h AB where g AB is the background metric, then under 
a gauge transformation x A — > x A + £ A the metric perturbation transforms as 

h A B -> h AB - g A cd B £ C - 9Bcd A t, C - £ C d c gAB- (97) 

Since a five- vector £4 has three scalar degrees of freedom £*, £ y and = Vi£ s , only four 
of the seven functions ($, T, W 7 , a, (3, x) are physical. This immediately tells us that we 
expect to be able to define four gauge invariant variables constructed from the metric alone. 
Let A denote ^ and A' denote Under a gauge transformation each of the variables 
transforms as 



n n 
t n n 

w^>w-z n + t 2 i\ 

b 2 ■ 

n z 

e^e 

b b' 

(98) 

It is then relatively easy to construct the following gauge invariant quantities 

®inv = $-a-a 13—, 

n n 

Y mv = T + (3' + -a + a- + (3-, 
t n n 

W mv = W -a' + t 2 k 

Vinv = ^ + -a + -(3, 
b b 

where a = a — -^x and (3 — (3 + -^rpx' ■ We then see that there is a special gauge defined 
byx = a = ( 5 = 0in which 

&inv = 

r = r 

wi„„ = w, 
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We define this to be five-dimensional longitudinal gauge and so we see that the gauge invari- 
ant variables equal the values of the metric perturbations in longitudinal gauge, in perfect 
analogy with four- dimensional cosmological perturbation theory. This gauge is characterized 
by being spatially isotropic in the x l co-ordinates but in general there will be a non-zero 
t — y component of the metric. 



Position of branes 

In general, the locations in y of the perturbed branes will be different in different gauges, 
and it is very important to understand this location in each case. Remarkably, in the case 
where the the brane matter has no anisotropic stress this is easy to establish. Start in the 
gauge a = x — 0. From the above transformation rules we can see that we can always go 
to this gauge using only £* and £ s transformation. This then leaves us with the freedom to 
perform any £ y transformation such that the position of each brane remains unperturbed. 
Then working out the Israel matching conditions we find that (3 on the branes is related 
to the anisotropic part of the brane's stress energy. So if we are considering only perfect 
fluids, for which the shear vanishes, then the Israel matching condition gives (3(y = y±) = 0. 
We can then go to longitudinal gauge (a — f3 — x — 0) with the transformation £ y — f3 
alone. But since /3 vanishes on the branes, so does £ y implying that the brane trajectories 
are unperturbed. So we see that for the special case of matter with no anisotropic stress 
the locations of the branes in longitudinal gauge are their unperturbed values y — y±. We 
can then infer the position of the branes in an arbitrary gauge by means of the above gauge 
transformations to be 

y = y±-P, (99) 

where y± are the background values. In particular in the class of Milne gauges we have 
defined in (57) the branes are located at 

y = y±- ^x'- (ioo) 



APPENDIX 3: BIRKHOFF'S THEOREM AND THE BACKGROUND METRIC 

The bulk geometry considered in this paper solves the five-dimensional Einstein's equa- 
tions sourced by a pure negative cosmological constant. For the background solution we 
restrict to solutions possessing cosmological symmetry on three dimensional spatial slices. 
In close analogy to the familiar situation for spherical symmetry in 3 + 1 dimensions, a 
Birkhoff-type theorem guarantees that in our case that away from the branes, the back- 
ground must take the form of either Anti-de Sitter (AdS) space-time, Schwarshild-AdS or 
AdS with a naked singularity. In each case the metric may be written as 

2 2 

ds 2 = (^- + k- ^Y l dr 2 - + k — ^)dT 2 + r 2 1%1 dx l dx\ (101) 
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where fi is the mass of the black hole, 7^ is the canonical metric on S* 3 , H 3 or E 3 , with 
k the corresponding spatial curvature, and L is the AdS radius defined by A = — 6Mf/L 2 
with M 5 the five-dimensional Planck mass. We are most interested in the case k = 0, for 
which it is useful to change variables from r to Y obtained by setting the first term in (101) 
to equal dY 2 , obtaining 

ds 2 = dY 2 - N(Y) 2 dT 2 + A(Y) 2 dx 2 , (102) 

where for AdS 

A{Yf = N(Y) 2 = exp[2y/L], (103) 
for Schwarzschild-AdS with a horizon at Y = 

A{Y) 2 = cosh(2F/L) and N(Y) 2 = gg^ , (104) 
and for AdS with a naked singularity at Y = 

MY? = ^(2YIL), N{Y? = ^1^. (1 05) 

For any configuration of branes possessing cosmological symmetry, even if the branes 
move the Birkhoff theorem guarantees that the bulk geometry takes one of the three forms 
above. 37,40 In our case, where the branes are ^-symmetric and have their tensions tuned 
to allow static empty brane solutions, the only bulk solution that is consistent with mov- 
ing branes is the Schwarzshild-AdS solution. Consequently this is the background five- 
dimensional metric we use in this paper. 

Technically, in order to study the perturbations it is much simpler if one changes coor- 
dinates to those in which the branes are static and the bulk is time-dependent. That it is 
always possible to choose such a coordinate system may be seen as follows. Start with the 
Birkhoff-frame metric (102) with A and N given by (104). First, change variables from Y 
to Z defined by dZ = dY/N, with Z chosen to be zero at the collision event, so that 

ds 2 = N 2 (-dT 2 + dZ 2 ) + A 2 dx 2 , (106) 

where N and A are now functions of Z. Defining lightcone co-ordinates T± = T± Z we have 

ds 2 = N 2 (—dT + dTJ) + A 2 dx 2 . (107) 

We now recognize that the form of this metric is invariant under the light-cone coordinate 
transformation, r ± y — f±(T ± Z), which takes the metric to the form 

N 2 

ds 2 = -jrjri-dr 2 + dy 2 ) + A 2 dx 2 . (108) 
/+/- 

Now we set t = ±e ±T , to describe the post- or pre-collision space-times respectively, and 
define t 2 n 2 (t, y) = N 2 /(f' + f'_) and b 2 (t, y) = A 2 to obtain 
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ds 2 = n 2 (t, y)(-dt 2 + t 2 dy 2 ) + b 2 {t, y)dx 2 , 



(109) 



which is the form used in this paper. 

We now show that we can always choose the functions f± to make the branes static in 
the new coordinates. To see this note that the new spatial coordinate 

y(T, Z)= l - (f + (T + Z)-f_(T- Z)) (110) 

itself satisfies the massless field equation in two dimensions. If the two brane trajectories 
are Z = Z±(T) in the T, Z coordinates, then it follows from the general theory of the 
wave equation that we can always solve (110) for arbitrary chosen y(T, Z) on two specified 
timelike curves Z = Z±(T). In particular we are free to choose constant values y = y + on 
the positive tension brane and y = y_ on the negative tension brane. Even after this choice 
there is additional coordinate freedom, since to determine the solution for y(T, Z) we need 
to specify additional Cauchy data, for example on a T =constant surface. 

In practice we find it is straightforward to solve these equations as a power series in t. 
The Israel matching conditions on the two branes in Birkhoff coordinates read 

tanh(2F ± /L) = (1 ± E^L) ^ _ N'\Y ± ){dY ± / dT)\ (111) 

where p± are the densities of matter or radiation on the branes. In our case, when only 
radiation is present, and we normalize the brane scale factors to be unity at collision (cf. 
Section III.C), we have p± = r±j A 4 (Y±). Equation (111) is a first order differential equation 
for the brane trajectories Y±(T), allowing them to be straightforwardly determined as Taylor 
series in T. Likewise we may solve explicitly for Z, 

^) = §(tan- 1 (x) + iln(^)), (112) 

where x 2 = cosh(2y/L), and hence obtain Z±(T) as a Taylor series in T. From (110) we 
obtain 

V± = I (U(T + Z ± (T)) -f.(T- Z ± {T))) , (113) 

which we may differentiate with respect to T, noting that the y± are constant, to obtain 

f' + (T + Z ± (T))(1 + V ± (T)) =f'_(T- Z ± (T))(1 - V ± (T)), (114) 

where V±(T) = (dZ±(T)/dT) are the brane velocities. These two equations may be simul- 
taneously solved as a power series in T with the ansatz f'±(z) = z^ 1 + /° + f±z + 

They are both trivially satisfied at order T -1 . At each subsequent power T n ,n > 
one obtains two equations which fix the two constants /±. Finally, writing f±(z) = 
c± + In z + f±z + f±z 2 /2 + . . ., with c± constants, we can write the equation for y± and take 
the limit T — > on the right hand side to obtain 

y ± = i(c + - c_) + i In (i±^) = i(c + - c_) + Of, (115) 
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where 6± are the rapidities of the positive and negative tension branes in the Birkhoff frame. 
Likewise we obtain (for t > 0) 

r = ^(c+ + c_) + lnt. (116) 

Setting t = ±e ±T as we do for t > or t < respectively, and choosing y + = — y_ 
(i.e. the Lorentz frame in which the branes have equal and opposite speeds), then fixes 
c + = — c_ = — \{6+ + 6 B ) = —6 B . Now one may invert the equations r ± y = f±(T ± Z) to 
express T + Z as a Taylor series in te y for t > (or te~ y for t < 0) and similarly T — Z as a 
Taylor series in te~ y (or te^). For example, post-collision one obtains 

T + Z = te y e W + 0(t 2 ), T - Z = te^e^ + 0(t 2 ), (117) 

equations which will be useful in Appendix 4. Hence we completely determine the metric 
functions n 2 and b 2 as Taylor series in te y and te~ y . Finally, by rescaling t and x we can 
also ensure that in the new coordinates, n(t,y) = 1 + 0(t) and b(t, y) = 1 + 0(t). 

As a check of this procedure, or indeed an alternative to it, one can directly solve Ein- 
stein's equations in the frame in which the branes are static. The extrinsic curvature is given 
by 

K^dx"dx u = ^-d y g^dx»dx v = ^-(-(n 2 )'dt 2 + (b 2 )'dx 2 ), (118) 
and so the Israel matching conditions 

1 ,„ 1 



^ = 2M| (T ^~3^ Ta) ' (U9) 
tell us that 



n' 



„n I T 3" ±T 2 p - (120) 

For the purposes of our analysis it will be convenient to define the Lorentz frame we work 
in to be that in which the y coordinates of the branes (their rapidities) are y± = ±y /2. 
Recall, we also define the parameters r± to be the the densities of radiation on each brane 
p± at collision, and we treat these as free parameters. Through a direct series solution of 
the five-dimensional Einstein equations, imposing the Israel matching conditions (121) at 
each order in t, we obtain the following solution for the background geometry near t — 0: 



b(t, y) = 1 + (bi sinhy + b 2 coshy)t + (e + e\ sinh2y + e 2 cosh2y)t 2 /2, 
n(t,y) = 1 + (disinhy + d 2 coshy)t + (k + h sinh 2y + k 2 cosh 2y)t 2 /2, (122) 



where the constant parameters are given by 
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fc = (12 + ^- r -»U (W2). 
&2 = + r -) cosech (V 2 )> 

4 = (4 - L2 ^- r - )) 3 ech(W2). 
d 2 = ~^( r + + r_)cosech(y /2), 

60 = 3^ (( " 6 + LV " )2 + (6 + LV+)2 + 2( " 6 + LV " )(6 + LV+) COSh2/0 ' 
+36(cosh2yo — l))(cosech|/o) 2 , 

ei = Y^(- 4 - £ 2 ( r + - r_))(r+ + r_)cosech(y ), 

e 2 = -^2"(( 24 + 4L2 ( r - - r +) + 2L 4 r + r_ + (-24 - 4L 2 (r_ - r+), 

+L 4 (ri + r 4 )) coshyo)(cosech?/o) 2 , 

A; = ;r^((21 + L 4 (r 2 + r 2 ) + 2(-12 + L 4 r_r + ) coshy + 3 cosh2y )(coseclM/ ) 2 , 
6L Z 

fci = -y^( r + + r_)cosechyo(5L 2 (r_ - r+) + 12sechy ), 

A; 2 = _i_(_24 + 10L 4 r + r_ + (24 + 5L 4 (r 2 + r 2 )) cosht/ )(cosech?/o) 2 . (123) 



APPENDIX 4: MEANING OF THE CONSTANTS d AND C 2 

The two arbitrary gauge constants C\ and c 2 in Eq. (67) parameterizing the violation of the 
asymptotically Neumann boundary condition (64) have a simple geometrical interpretation: 
they describe the displacement of the collision event in the T, Y plane. Recall that in 
the Neumann gauge, discussed in Section V.E, the brane trajectories are unperturbed and 
are are described by the equations y = y± = constant. If we now gauge transform to an 
asymptotically Neumann gauge, in which the normal derivatives 7i -1 £ -1 x'(£, y±) deviate from 
zero at order t as in (67), we see that the gauge transformation from conformal Newtonian 
gauge to the Milne gauge we are in involves a divergent y coordinate displacement of £ y = 
—q 2 x'/t 2 , which tends to — (cisinhy + c 2 coshy)/t plus a finite part as t tends to zero. 
If c\ = c 2 = 0, then the perturbation in the brane y coordinates, £ y (y±) is finite. The 
rapidities of the two branes are perturbed, but the collision event itself is still simultaneous 
as in Neumann gauge. 

In the remainder of this Appendix we provide a geometrical interpretation of the two 
constants c\ and c 2 , showing that they parameterize the displacement of the brane collision 
event away from its background location in the embedding coordinates T, Y, at each x. 

If we start from Neumann gauge with c\ = c 2 = 0, we may introduce c\ and c 2 via the 
following gauge transformation, 

£ s = [c\ cosh y + c 2 sinh y)t, 
£ y = — (ci sinh y + c 2 cosh y) , 

V 
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= (ci cosh y + c 2 sinh y) . 



This is part of the gauge freedom described by the solutions to eq.(59) and (60). Although 
£ y diverges near t — 0, this is merely a reflection of the singular nature of the Milne (£, y) 
coordinate system. In terms of the Birkhoff frame T, F coordinates defined in Appendix 3, 
we find 

dT dT 
ot oy 

5Y = f e + ?L e . 

ot ay 

Then using (117) given in Appendix 2 and (124) one infers the displacement of the collision 
event 

ST = (ci cosh W - c 2 sinhP") , 

SY = N{Y C ) (ci sinh^ - c 2 cosh^) , (124) 

independent of y and hence holding for both branes. Here N(Yc) is the value of the lapse 
function (given in Appendix C) at the collision value of Y in the Birkhoff frame, and 9 B is 
the mean rapidity of the two branes in that frame. Therefore all the gauge transformation 
(124) does is to move the collision event around by an arbitrary finite displacement in the 
T, Y plane. 
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